G 

(0 

O 


D DO 

7EJME- 

JUL  271919 

I \ VetEjlG 

r> 


Technical  Report  387 


©: 

00 

i> 

rH 

© 


'a 

s 


NAVAL  BLUE-GREEN  SINGLE-PULSE 
DOWNLINK  PROPAGATION  MODEL 

Technical  Advisor  to  the  Blue-Green 
Optical  Communication  Program 
Joint  Coordinating  Committee 

1 January  1979 

Final  Report:  January  1978  — December  1978 


Prepared  for 
Naval  Electronic  Systems  Command 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CALIFORNIA  921S2 

7 o n n' 


Q O 


^11 


NOSC  TR  387 


( 


NAVAL  OCEAN  SYSTEMS  CENTER.  SAN  OIEQO.  CA  921S2 


AN  ACTIVITY  OF  THE  NAVAL 

RR  GAVAZ2I , CAPT,  USN 

Commander 


MATERIAL  COMMANO 
HL  BLOOD 

Technical  Director 


ADMINISTRATIVE  INFORMATION 

The  analytical  model  for  blue-green  pulse  propagation  reported  within  this  document 
was  the  result  of  the  combined  efforts  of  Dr.  P.  J.  Titterton  and  H.  E..JS»figfley..of  GTE 
Svlvania  R D.  Anderson[  R,  F.  Howarth,  M.  E.  Hyde  and  L.  B.  Stotts  of  the  Naval  Ocean 
Systems  Center,  and  W.  R.  Stone  of  the  MEGATEK  Corporation.  The  authors  would  like 
to  thank  Mr.  M.  L.  Parker,  Jr.,  Naval  Electronics  System  Command,  Code  3102,  for  his 
initiation,  interest  and  support  of  this  work. 


Released  by 

LCDR  S.  E.  Dollar,  Head 
Communications  Research 
and  Technology  Division 


Under  Authority  of 
R.  H.  DuBois,  Head 
Communication  Systems 
and  Technology  Department 


4 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whan  Du*  tnlarag) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

NOSC  Technical  Report  387  (TR  387) 

2.  GOVT  ACCESSION  NO. 

1.  RECIPIENT'S  CATALOG  NUMBER 

/T7  

1 


TITLE  (end 


JiAVALBLUE-GREEN SJNGLEWLSE  DOWNLINK 
j^OPACATlObf MODEL  , ~ -T- 


"rr^T 

Final  Repast.  Jan, 


nil Report . Jan  4HB-Dec/M78  . 
NFOmllN^N^JtFoSTmiffDtN 


7.  A'jTHONf*; 

Technical  Advisor  to  the  Blue-Green  Optical  Communication 
Program  Joint  Coordinating  Committee 


t.  FfAFORMINO  ORGANIZATION  NAM*  ANO  ADDRESS 

Naval  Ocean  Systems  Center  7 

San  Diego,  CA  92152 


■ I.  CONTROLLING  OFFICE  NAMC  ANO  ADDRESS 

Naval  Electronic  Systems  Command 
Washington,  DC  20360 


T,  PROJECT,  TASK 
C WORK  UNIT  NUMBERS 


Ol 


MONITORING  AGENCY  NAME  A AODRCSSfU  dlllarant  Item  Controlling  Otllco)  IS.  SECURITY  CLASS,  (al  i’ll*  raporl) 

Unclassified 


17.  distribution  stateme 


is.  supplementary  no 


P.  J./Titt-e rt.cn, v"'  "•  ~ ./Sweeney 
R.  P.  'nrfcx’con,  J. . F . /Rower  th 
~./!!ydc,  , - - — — — - — ■ 

• 


If.  KEY  WO  ADS  (Continue  on  «vfri«  aide  If  ne cm  a aery  end  Identify  by  block  number) 


JOX^AIITHACT  (Cenllnue  on  tewetee  cldc  II  nacaaaery  and  Identify  by  block  number) 

-^-This  technical  report  describes  the  official  Navy  model  for  signal  and  noise  characteristics  of  a single  blue-greer 
optical  pulse  propagated  from  a satellite  laser  source  to  a submerged  receiver  in  all  weather  conditions. 

A mathematical  model  is  developed  in  two  group?:  the  signal  group  relates  the  instantaneous  received  optical 
signal  power,  while,  the  noise  group  relates  the  total  ,ione  equivalent  optical  power  to  the  noise  source.  Both 
groups  of  the  model  account  for  propagation  path  and  receiver  characteristics;  and,  they  allow  calculation  of 
the  signal-to-noise  ratio  for  any  blue-green  single-pulse  downlink  from  a satellite  laser  source  to  a submerged 
receiver.  The  relative  technical  uncertainty  in  accuracy  of  both  equation  groups  has  also  been  included. 


00  , 1473 


COITION  of  I NOV  ai  is  obsolete 

S/>4  0102-LF.014.e60! 


UNCLASSIFIED  \ 

SECURITY  CLASSIFICATION  O*  This  RAaE  (Whtn  Duo  BMata*) 

~39  3 


i 

i 


t 


SUMMARY 

OBJECTIVE 

Develop  a mathematical  model  for  the  noise  and  signal  characteristics  involved  in  the 
propagation  of  a single  blue-green  optical  pulse  from  a satellite  laser  source  to  a submerged 
receiver  under  all  weather  conditions. 

RESULTS 

A mathematical  model  is  developed  in  two  groups:  the  signal  group  relates  the 
instantaneous  received  optical  signal  power,  while,  the  noise  group  relates  the  total  noise 
equivalent  optical  power  to  the  noise  source.  Both  groups  of  the  model  account  for 
propagation  path  and  receiver  characteristics;  and,  they  allow  calculation  of  the 
signal-to-noise  ratio  for  any  blue-green  downlink  from  a satellite  laser  source  to  a 
submerged  receiver.  The  relative  technical  uncertainty  in  accuracy  of  both  equation 
groups  has  also  been  included. 
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1.0  INTRODUCTION 


The  Navy  has  long  recognized  that  the  so-called  "blue-green  window”  for  optical 
propagation  through  sea  water  offers  a potential  communication  channel  for  information 
transfer  to  submarines  at  operational  depth  and  speed.  This  recognition  has  spurred  contin- 
uing efforts  to  develop  a technology  base  in  blue-green  propagation  and  devices.  This 
report  is  the  initial  documentation  of  the  official  Navy  model  for  the  signal  and  noise 
characteristics  of  a single  blue-green  optical  pulse,  propagating  from  a satellite  laser  source 
to  a submarine  receiver  in  all  weather  conditions. 

The  model  is  divided  into  two  groups  of  equations:  the  signal  group  relates  the 
instantaneous  received  optical  signal  power,  while,  the  noise  group  relates  the  total  noise 
equivalent  optical  power  to  the  source.  Both  groups  of  the  model  account  for  propagation 
path  and  receiver  characteristics:  and,  they  allow  calculation  of  the  SNR  for  any  single-pulse 
blue-green  downlink  from  a satellite  laser  source  to  a submerged  receiver. 

Figure  I shows  all  the  input  parameters  for  each  portion  of  the  signal  group.  The 
laser  pulse  originates  at  the  satellite,  propagates  through  the  atmosphere  (including  whatever 
clouds  are  present),  the  air-sea  interface,  the  ocean  water,  and  is  finally  detected  by  the  sub- 
merged receiver.  With  noise,  some  of  the  background  light  sources  (sun  and  moon)  follow  a 
like  path,  while  others  (blue  skylight  and  star/zodiacal  light)  are  not  point-like  beam  sources 
and  must  be  treated  differently.  Finally,  the  bioluminescent  light  originates  in  the  water 
itself  and  is  not  effected  by  atmospheric  properties. 

In  this  report,  sections  2.0  and  3.0  relate  the  model  and  are  organized  in  the  follow- 
ing manner.  First,  the  method  of  approach  used  in  each  group  is  briefly  discussed  and  a 
detailed  flow  chart  (figs  3,  20  respectively)  showing  the  interrelations  between  all  the 
equations  is  illustrated.  The  second  subsections  (2.2  and  3.2)  define  and  discuss  all  input 
information  needed  to  perform  the  calculations.  The  third  and  final  subsections  (2.3  and 
3.3)  contain  the  derivations  and  justifications  of  all  equations  used.  These  equations  and 
derivations  are  discretely  organized  in  a modular  fashion,  so  future  revisions  (e.g.,  cloud  and 
water)  may  be  made  without  major  dislocations  to  the  model. 
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Figure  1.  Parameters  for  the  single-pulse  downlink  model. 

2.0  SIGNAL  PROPAGATION  GROUP 


This  section  discusses  the  equation  group  for  the  optical  propagation  of  a single 
blue-green  signal  pulse  front  a satellite  to  a submerged  receiver.  The  output  of  the  signal 
group  is  combined  with  that  of  the  noise  group,  to  define  the  signal  to  noise  ratio  (SNR) 
for  a given  communication. 

2.1  OBJECTIVE,  STRUCTURE  AND  FUNCTION 

This  section  describes  the  basic  approach  used  with  the  detailed  single  pulse 
downlink  signal  propagation  group  and  presents  a flow  chart  showing  the  interrelationships 
of  the  sub-models  and  their  required  inputs. 

2.1.1  Philosophy  of  Approach 

The  signal  group  is  organized  in  a modular  fashion,  so  that  the  effect  of  each  portion 
of  the  path  is  evident.  In  addition,  as  further  experiments  and  analyses  are  undertaken, 
pieces  of  the  model  may  be  upgraded  without  requiring  extensive  modification  of  the  total 
model.  The  signal  calculation  group  considers  the  following  three  properties  of  the  signal 
and  separately  models  the  effects  of  the  propagation  path  on  each: 

(a)  the  energy  transmission  from  the  satellite  transmitter  to  the  submerged 
receiver, 

(b)  the  distribution  of  the  signal  radiance  at  the  submarine  receiver,  and 

(c)  the  pulse  shape  and  width  at  the  submarine  receiver. 

These  three  properties  are  then  combined  to  yield  the  instantaneous  received  optical  power 
at  the  surface  of  the  submerged  photodetector. 

This  approach  does  not  attempt  to  treat  all  possible  cloud  conditions.  Rather,  a 
break  point  is  established  at  a minimum  optical  thickness  of  10.  Below  that  value,  one  set  of 
sub-models  is  assumed  to  apply,  while  above  it  a different  set  applies.  In  many  cases,  these 
sub  models  do  not  correspond  at  rQpt  s optical  thickness  - 10,  and  so  the  overall  model 
should  only  be  used  for  rQpt  <10  and  rQpt  > 10.  (Future  analysis  and  experiments  on  the 
“multiple-forward  scattering”  region  should  enable  the  sub-models  to  be  upgraded,  remov- 
ing this  inconsistency. 

Finally,  simple  analytic  forms  are  assumed  for  at-present  unknown  functions  such  as 
the  radiance  distribution,  and  pulse  shape  and  width.  This  enables  production  of  analytic 
results  (except  for  the  receiver  axis  off-set  from  the  beam  axis  of  the  incident  radiance), 
which  are  an  aid  to  a physical  understanding  of  the  overall  propagation  problem. 

2. 1 .2  Signal  Group  Structure 

A schematic  of  the  overall  signal  propagation  group  is  presented  in  figure  2.  Given 
the  input  parameters,  the  pulse  width,  shape,  and  angular  and  radiance  distribution  are 
derived.  Then,  using  the  r igular  and  radiance  distribution,  and  the  input  data,  the  energy 
transmission  of  the  path  is  calculated.  Finally,  given  the  energy  transmission  and  the 
received  pulse  shape,  the  signal  (instantaneous  received  power)  is  calculated. 


Figure  2.  Schematic  of  single-pulse  downlink  signal  propagation  group. 

Figure  3 is  a detailed  flow  diagram  showing  the  calculations  that  must  be  performed 
to  arrive  at  the  required  output.  The  following  five  steps  are  used  to  exercise  the  signal 
propagation  group: 

• The  input  parameters  are  listed  in  the  seven  ellipses  on  the  left  hand  side  of  the 
figure,  including  source,  clear  atmosphere,  cloud,  cloud  to  water,  air-water  interface, 
water,  and  receiver  parameters.  (The  symbols  are  defined  in  the  glossary  and  also  in 
the  input  discussion  in  section  2.2). 

• The  primary  calculation  equations  are  represented  by  the  rectangular  boxes.  Within 
each  box  is  the  symbol  for  the  parameter  to  be  calculated  and  the  equation  number 
(from  section  2.3)  for  the  equation  to  be  used  to  calculate  the  parameter.  The  first 
quantity  calculated  is  the  cloud  optical  thickness,  tq pt  , since  this  determines  the 
equations  to  be  used  to  calculate  many  other  parameters.  Whenever//  appears  in  a 
rectangular  box,  the  equation  number  preceding  it  refers  to  r0pt  > 10,  while  the 
equation  number  following  it  refers  to  r0pt  < 10-  Hence,  given  the  value  of 
r0pt  the  rest  of  the  models  to  be  used  are  specifically  determined. 

• The  second  set  of  calculations  performed  are  of  three  types: 

(a)  path  transmission,  including  ra,  rc,  rcw,  raw  and  rw; 

(b)  pulse  width  and  shape,  including  Atc,  Atcw,  Atw  and  t^j. 

(c)  angular  and  radiance  distribution,  including  f^,  f^yy,  fyy,  <PQ  and  f(0^  <pQ,  6), 

• The  path  transmission,  angular  and  radiance  distribution,  source  and  receiver  parameters 
are  then  used  to  calculate  the  received  energy,  Er. 

• The  received  energy,  pulse  width  and  shape  are  used  to  calculate  the  instantaneous 
received  power  Pj^(t). 

The  computer  program  to  be  developed  will  follow  the  flow  chart  in  figure  3. 
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INPUT 

PARAMETERS 


PB(t)  : (S-33) 


(NOISE) 

NEP,-, 


LEOENO: 

PARAMETER: 

THICK  CLOUD  E'Q'N// 

THIN  CLOUD  E'Q'N 

PARAMETERS  DEFINED  IN  GLOSSARY 
INPUT  DISCUSSION  IN  SECTION  2.2 
EQUATIONS  DEVELOPED  IN  SECTION  23 


Figure  3.  Flow  diagram  of  single-pul  ie  downlink  tignal  propagation  group. 
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2.2  INPUT  FOR  SIGNAL  CALCULATION 


This  section  discusses  the  form  of  the  required  inputs  to  the  single  pulse  downlink 
signal  propagation  model,  in  terms  of  the  seven  categories:  source,  clear  atmosphere, 
cloud,  cloud  to  water,  air-water  interface,  water  and  receiver. 

2.2.1  Source 


The  required  source  parameters  are: 


Symbol 

q 


“P 

yT 

R 

0T 


Description 

A parameter  describing  the  ability  of  the  satellite  transmitter  to 
correct  for  the  geometric  zenith  angle  spreading  of  the  spot,  q * 0 
implies  the  spot  remains  the  same  area,  independent  of  zenith 
angle,  while  q = 1 means  the  spot  grows  naturally  with  zenith  angle. 
Hence,  0 < q < 1 . 

Energy  per  pulse  of  the  transmitter  laser. 

Transmission  of  the  transmitter  optical  chain. 

Range  from  the  satellite  transmitter  to  the  submarine  receiver. 


Units 


Full  angle  beam  divergence  to  the  e 
transmitter  beam. 


-2 


irradiance  points  of  the 


Joules 


metres 

radians 


2.2.2  Clear  Atmosphere 

The  required  parameters  are: 

Symbol  Description 

b Effective  clear  atmosphere  optical  thickness  such  that  for  a 70 

percent-zenith  transmission,  b = 0.357 

0-  In-air  transmitter  zenith  angle 


Units 


radians 


2.2.3  Cloud 


The  required  parameters  are: 

Symbol  Description 

T Geometric  or  physical  thickness  of  the  cloud 
ac  Average  extinction  coefficient  of  the  cloud 

q A parameter  describing  the  ability  of  the  satellite  transmitter  to 

correct  for  the  geometric  zenith  angle  spreading  of  the  spot, 
q = 0 implies  the  spot  remains  the  same  area  independent  of 
zenith  angle,  while  q = 1 means  the  spot  grows  naturally  with 
zenith  angle.  Hence  o < q < 1 . 

<cos  8>  The  average  value  of  the  cosine  of  the  scattering  angle  for 
single  scattering  within  the  cloud 

<t> <.  In-air  transmitter  zenith  angle 


Units 

metres 

(metres)-' 


radians 
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Symbol 

Description 

Units 

(jJ- 

o 

The  single  scattering  albedo  of  the  cloud  particles,  or,  the  ratio 
of  the  probability  of  scattering  to  the  probability  of  extinction 
for  a single  scattering  event. 

6 

The  average  angle  for  single  scattering  within  the  cloud 

radians 

2.2.4  Cloud  to  Water 

The  required  parameters  are: 

Symbol 

Description 

Units 

H 

Distance  from  cloud  base  to  water  to  surface 

metres 

R 

Range  from  the  satellite  transmitter  to  the  submarine  receiver 

metres 

Full  angle  beam  divergence  to  the  e“^  irradiance  points  of 
transmitted  beam 

radians 

2.2.5  Air-Water  Interface 

The  required  parameters  are: 

Symbol 

Description 

Units 

*s 

In-air  transmitter  zenith  angle 

radians 

V 

Surface  wind  speed 

metres/ 

second 

n 

Water  index  of  refraction 

2.2.6  Water 

The  required  parameters  are: 

Symbol 

Description 

Units 

ki 

Diffuse  attenuation  coefficient  of  the  i’th  water  layer 

(metres)- 

Di 

Thickness  of  the  i’th  water  layer 

metres 

In-water  transmitter  zenith  angle 

radians 

*si 

Root-mean-square  angle  for  a single  scattering  event  in  the 
water 

radians 

s 

Scattering  coefficient  for  the  entire  water  path 

(metres)- 

n 

Water  index  of  refraction 

D 

Depth  of  the  submarine  receiver 

metres 

2.2.7  Receiver 

The  required  parameters  are: 
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Symbol 

Description 

Units 

0R 

Half-angle  of  the  receiver  field  of  view 

radians 

5 

Off-set  angle  between  the  in-water  signal  beam  axis  and  the 
receiver  optical  axis 

radians 

D 

Depth  of  the  submarine  receiver 

metres 

tr 

Transmits  on  of  the  receiver  optical  chain 

d 

Diameter  of  the  receiver  optical  aperture 

metres 

2.3  SIGNAL  CALCULATION  MODULES 

This  section  develops  all  the  equations  used  in  the  calculation  of  the  instantaneous 
received  signal  power.  Sections  2.3.1,  2.3.2,  2.3.3,  2.3.4,  2.3.6  and  2.3.7  consider  the  path 
transmission  of  the  energy.  Sections  2.3.5  and  2.3.8  consider  the  angular  effects  and  the 
distribution  of  the  received  radiance.  Section  2.3.9  considers  the  received  pulse  shape  and 
width.  In  each  of  these  sections,  after  the  equations  are  developed  they  are  evaluated  for 
typical  cases  in  both  tables  and  figures.  Section  2.3.10  combines  the  previous  results  to 
obtain  the  received  energy  and  the  optical  signal  power. 


2.3. 1 Clear  Atmospheric  Transmission 


In  the  absence  of  any  clouds  or  aersols,  the  clear  atmospheric  transmission  is 
described  by  the  term  ra.  Using  the  approximate  AFCRL  model  (ref  1),  the  signal  zenith 
angle  dependence  is  given  by: 


ra  = exp  [-bsec  0S]  , 

for 


(S-l) 


ra  = signal  clear  atmospheric  transmission, 

b = effective  clear  atmospheric  optical  thickness,  and 
0S  = signal  zenith  angle. 

The  typical  value  of  b is  determined  from: 
ra(0s  = O)  = 0.7  = exp(-b)  , 


ot 


b = 0.357  . 

Table  1 and  figure  4 show  the  values  of  ra  as  a function  of  signal  zenith  angle. 


I . AFCRL-7 1 -0279,  Optical  Properties  of  the  Atmosphere  (Revised),  by  RA  McClatchey,  RW  Fenn, 
JEA  Selby,  FE  Volz,  and  JS  Garing,  10  May  1971 . 
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Table  1.  Typical  clean 
0S,  Single  Zenith  Angle  (degrees) 


* oa 

z 

o 

8 0.7 

S 

Z 0.8 

e 

o os 
c 


i °-3 


10  IS  20  25  30  35  40  45  50  55  60  65  70  75  80  85 


SIGNAL  ZENITH  ANGLE,  <Pt 

Figure  4.  Typical  clear  atmospheric  tranimiMion  (b  = 0.357), 

2.3.2  Goud  Energy  Transmission 

Insofar  as  they  affect  optical  propagation,  clouds  may  be  categorized  as  negligible 
to  very  thin,  thin,  and  thick.  There  are  no  verified  experimental  results  that  treat  any  one 
of  these  three  classes.  Since  we  are  discussing  an  energy  transmission  here,  and  most  analyses 
and  partial  experiments  are  in  terms  of  the  transmitted  irradiance  (energy/second/area), 
there  are  few  analytic  expressions  for  the  cloud  energy  transmission. 

Our  decision  is  to  treat  only  two  classes,  thick  clouds  and  nearly  clear  weather. 

We  do  this  with  the  understanding  that  evaluating  the  model  near  the  transition  point 
will  yield  results  that  are  incorrect  in  principle,  and  so  only  for  small  and  large  optical 
thickness  (defined  below)  should  the  overall  model  be  expected  to  apply. 
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We  propose  to  adopt  the  multiple-scattering  Monte-Carlo  derived  diffusion-like 
expression  of  Bucher  (ref  2)  and  Van  der  Hulst  (ref  3)  for  thick  clouds,  so  that  at  zenith: 

1.69 (S-2) 

T°  TopX(\-<coS6»+lA2 

where 

rc  = signal  energy  transmission  through  the  cloud, 

ropt  = optical  thickness  of  the  cloud,  and 
<cos  6>  - mean  cosine  of  the  scattering  angle. 

In  discussing  this  equation,  Bucher  has  estimated  that  the  numerical  parameters 
( 1 .69  and  l .42)  are  uncertain  to  within  factors  of  2X. 

The  optical  thickness  of  the  cloud,  for  a homogeneous  cloud,  is  given  by: 

ropt  = Toc  > (S‘I * 3) 

where 


T = geometrical  thickness  of  the  cloud,  and 
ac  = mean  extinction  coefficient  of  the  cloud. 

As  a tvpical  example,  for  a strato-cumulus  cloud,  we  might  have  T = 1200  m and 
oc  * 0.04  m , so  r0pt  = 48.  For  such  a dense  cloud,  <cos  0>  = 0.83  (0  * 34°);  and  so, 

t = — =0.18  . 

c 48(1  -0.83+1.42 

The  zenith  cloud  energy  transmission,  as  a function  of  optical  thickness,  is  presented 
in  table  2 and  figure  5 for  <cos  9>  = 0.83. 

From  the  form  of  equation  (S-2)  it  cannot  apply  ropt  = 0,  and  the  lower  limit  of 
optical  thicknesses  at  which  it  does  apply  is  still  to  be  determined.  A linear  fit  is  provision- 
ally adopted  for  ropt  < 1 0,  so  that 

rc  = 1 - 0.046  ropt  , Topt  <10 

<cos  6>  = 0.83 

This  equation  is  evaluated  in  table  3,  and  figure  5. 


I 

2.  EA  Bucher,  Computer  Simulation  of  Light  Pulse  Propagation  for  Communication  Through  Thick 
Clouds,  Appl  Optics,  Vol  12  (10),  p 2391 -2400,  October  1973. 

3.  RE  Danielson,  DR  Moore,  and  HC  Van  de  Hulst,  The  Transfer  of  Visible  Radiation  Through  Clouds, 

j J Atmos  Sci,  Vol  26  (9),  p 1078-1087,  September  1969. 

i 
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CLOUD  ENERGY  TRANSMISSION, T 


Table  2.  Typical  “thick"  cloud  zenith  signal  energy 
transmission  (<cos  0>  * 0,83). 


r0pt  , Optical  Thickness 

rc,  Energy  Transmission 

10 

0.54 

20 

0.35 

30 

0.26 

40 

0.21 

50 

0.17 

60 

0.15 

70 

0.13 

80 

0.11 

90 

0.10 

100 

0.09 

Figures.  Thin  and  thick  cloud  energy  trammUsion  versus  optical 
thickness  for  <cos  d>  = 0.83. 
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Table  3.  “Thin”  cloud  zenith  signal  energy  transmission  (matched 
to  the  thick  cloud  expression  at  rQpt  = 10  for  <cos  6>  = 0.83). 


ropt’  Optica]  Thickness 

rc,  Energy  Transmission 

0 

1 

2 

0.91 

4 

0.82 

6 

0.72 

8 

0.63 

10 

0.54 

The  zenith  angle  dependence  for  thick  clouds  has  also  been  modeled  by  Bucher  and 
Van  der  Hulst  (refs  2,3).  It  also  has  no  experimental  verification.  From  their  tabuiar  results, 
we  find 

rc  a M°>  A(*s>  ’ 

when 

A(<*»s)  = 1 .69  - 0.55 1 3 + 2. 7 1 73  <j>l  - 6.9866  <p f 

+ 7.1446^-3.4249^  + 0.6155^  , (S-4) 

neglecting  the  effects  of  the  “spreading  out”  of  the  beam  energy  with  zenith  angle.  Since 
this  “spreading  out”  may  be  corrected  at  the  satellite  terminal,  we  model  it  as  a transmission 
factor, 

(cos  <t>s)q  , 0 < q « 1 ~~  . 


The  complete  thick  cloud  signal  energy  transmission  (including  the  effects  of  the 
single  scatter  albedo  (ref  4),  cjq.  % 1)  is  given  by: 


4.  The  effect  of  coq  ^ 1 is  taken  from  Appendix  B,  equation  14,  of  S.  Karp,  “A  Test  Plan  for 
Determining  the  Feasibility  of  Optical  Satellite  Communications  Through  Clouds  at  Visible 
Frequencies,”  NOSC  TN  279.  July  1 , 1978. 
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and  0 < q < 1,  with  the  exact  value  of  q dependent  on  the  satellite  optical  design. 


For  thin  clouds  the  diffusion,  like  A(^s)  dependence,  ought  not  to  apply.  Therefore, 
in  the  < 10  regime  we  assume  a full  cosine  dependence  due  to  the  extra  path  length 
in  the  cloud ; or, 


1 - 0.885rOpt 

r 169  fee 

>s*l  q+1 

rc 

110(1  -<cos$»  + 1.42  J 1 L" 

)S  <PSJ 

forTopt  < 10  • 

For 

ropt 

= 0,  we  use  rc  = J 

cos  0SJ  9 . 

(S-5b) 

For  <ps  % 0,  these  models  are  discontinuous  at  = 1 0.  The  difference  in  value  is 
about  a factor  of  0.58  at  #s  = 70  degrees.  At  this  stage  in  our  knowledge  of  cloud  propaga- 
tion, we  do  not  feel  that  such  a factor  is  of  prime  importance  and  shall  ignore  this  discrep- 
ancy until  an  experimentally  verified  model  replaces  it. 

Table  4 and  figure  6 show  the  zenith  angle  dependence  of  the  signal  energy  trans- 
mission for  both  ropt  regimes,  with  q = 0;  i.e.,  satellite  optics  fully  compensating  for  the 
zenith  angle  beam  irradiance  spread. 

Table  4.  Zenith  angle  vs  signal  energy  transmission  (normalized  to  = 0) 


0$,  Signal  Zenith  Angle 

Thick  Cloud  Dependence 
(Topt  > '0) 

Thin  Cloud  Dependence 
('opt  < 1°) 

0 

1 

1 

10 

0.97 

0.98 

20 

0.96 

0.94 

30 

0.92 

0.87 

40 

0.86 

0.77 

50 

0.78 

0.64 

60 

0.69 

0.5 

70 

0.58 

0.34 

80 

0.44 

0.17 

H 
i i 

[i 

1 1 


i 


* \ 
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Figure  6.  Thin  and  thick  cloud  zenith  angle  dependence  of 
cloud  transmissions  normalized  to  zenith. 


2.3.3  Cloud  to  Water  Energy  Transmission 

The  energy  transmission  of  an  assumed  clear  propagation  path  from  the  bottom  of 
the  cloud  to  the  water  surface  is  expressed  as  For  large  diameter  emerging  spots  and 
representative  clouds  the  energy  transmission  of  this  part  of  the  propagation  path  is  very  high. 

We  consider  an  area  element  on  the  cloud  bottom  given  by  r dr  d <p,  assuming  circular 
symmetry  and  that  our  observation  point  is  well  away  from  the  beam  edge. 

Then,  assuming  a thick  cloud,  the  energy  is  emerging  Lambationally,  and  taking 

H = distance  from  cloud  base  to  water. 


range  to  satellite, 


6j  = full  angle  e~“  irradiance  beamwidth  (<1°);  and 


the  beam  is  so  large  at  the  cloud  that  negligible  additional  in-cloud  induced  spreading  occurs. 
The  energy  transmission  is  given  by 


rCW  ■ 


-i  / 


> R«t  « 

dip  J r d r x j “tilt”  of  emitting  surface  areaj 
o o 

Range  from  surface  area  to  water  surface  “receiver”  | 

“tilt”  of  surface  "receiver’ 


21 


J 


or 


,2,  R0T/2 

rCW = " / d4>  J r d r 

* J o o 

H 


H 

1 

[h2ht2]* 

[h2  ♦ ,2]  * 

[H2  + r2]* 


thus, 


[H/(R  0t/2)1  2 

Tcw  - 1 - : ~ ’ r°Pl  > 10  • 1 


1 + [H/(R  dj/2  ] 


(S-6a) 


Equation  (S-6a)  is  evaluated  in  table  5 and  figure  7 as  a function  of  the  ratio  of  the 
cloud  base  height,  H,  to  the  emerging  spot  radius,  Rfl-p/2. 


Table  5.  Thick  cloud-to-water  surface  signal  energy  transmission. 


/R 6j  Cloud  Base  Height 

Energy  Transmission 

H / 2 Emerging  Spot  Radius 

0 

1 

0.05 

0.998 

0.1 

0.99 

0.15 

0.98 

0.2 

0.96 

0.25 

0.94 

03 

0.92 

0.35 

0.89 

0.4 

0.86 

0.45 

0.83 

0.5 

0.8 

IR6j 

For  typical  situations,  H j — y < 0.2  and  so  this  is  a very  small  effect  for  thick 

clouds. 

For  thin  clouds,  the  net  transmission  is  even  higher  for  large  incident  beams,  since 
the  energy  is  not  emerging  isotropically,  but  does  retain  its  emitted  angular  sense.  We 
therefore  take 


rCw=  !-  fort-opt  < 10  • 


(S-6b) 


We  further  assume  that  there  is  no  zenith  angle  dependence  for  tq^  in  either  class. 
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RATIO  OF  CLOUD  BASE  HEIGHT  TO  EMERGING  SPOT  RADIUS,  H/l I 

Figure  7.  Thick  cloud  to  water  surface  energy  transmission. 

2.3.4  Air-Water  Interface  Transmission 

The  signal  energy  transmission  of  the  air-water  interface  is  composed  of  two  factors; 

(S-7) 


Taw  (TawP  * ^raw2^’ 


where 


raw  = tota*  energy  transmission  of  air-water  interface,  with 

Tawl  = a*r*water  interface  transmission  due  to  index  of  refraction  discontinuity;  and 

Taw2  = air'water  interface  transmission  due  to  foam  and  streaks  on  the  sea  surface. 

This  section  treats  rawj  while  raw2  is  discussed  in  section  2.3.6. 

For  thin  clouds  and  clear  weather  (ropt  < 10),  we  assume  the  signal  beam  retains  its 
angular  sense.  Gordon  (ref  S)  has  estimated  tne  value  of  raw  j as  a function  of  surface  wind 
speed  and  signal  zenith  angle,  for  the  Cox  and  Munck  (ref  6)  wave  slope  distribution,  as 
shown  in  table  6 and  figure  8, 

For  thick  clouds  (rCpt  > 10),  the  energy  is  modelled  as  being  incident  on  the  sea 
surface  from  the  entire  hemisphere.  Using  the  Fresnel  refraction  equation,  and  neglecting 
wave  effects,  we  find 


5.  J Gordon,  Directional  Radiance  (Luminescence)  of  the  Sea  Surface,  S10,  ref  89-20,  October  1969. 

6.  C Cox  and  W Munck,  Statistics  of  the  Sea  Surface  Derived  from  Sun  Glitter, ) Marine  Res, 

Vol  13  2,  1954. 
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Fable  6.  rflW  Time-averaged  downlink  air-sea  interface  transmittance 
(for  thin  clouds,  rQpt  < 10). 


Signal 


v.  Wind  Speed 


Zenith 

Angle 

0 

1.03 

2.06 

4.12 

7.21 

10.3 

13.4 

16.5 

19.6  m/iec 

in  Air 

0 

2 

4 

8 

14 

20 

26 

32 

38  knots 

0 

0.979 

0.977 

0.976 

0.974 

0.970 

0.967 

0.963 

0.960 

0.9S6 

5 

0.975 

0.974 

0.972 

0.970 

0.966 

0.963 

0.959 

0.956 

0.952 

10 

0.964 

0.962 

0.961 

0.959 

0.955 

0.9S1 

0.948 

0.944 

0.941 

15 

0.945 

0.944 

0.943 

0.940 

0.936 

0.933 

0.929 

0.926 

0.922 

20 

0.920 

0.918 

0.917 

0.914 

0.910 

0.907 

0.903 

0.899 

0.896 

25 

0.887 

0.885 

0.884 

0.881 

0.877 

0.873 

0.870 

0.866 

0.863 

30 

0.847 

0.845 

0.844 

0.841 

0.837 

0.833 

0.829 

0.826 

0.822 

35 

0.800 

0.798 

0.797 

0.794 

0.790 

0.786 

0.782 

0.779 

0.775 

40 

0.747 

0.745 

0.743 

0.741 

0.736 

0.733 

0.729 

0.725 

0.722 

45 

0.687 

0.685 

0.684 

0.681 

0.677 

0.673 

0.669 

0.666 

0.663 

50 

0.620 

0.619 

0.617 

0.615 

0.611 

0.608 

0.605 

0.602 

0.599 

55 

0.548 

0.546 

0.545 

0.543 

0.540 

0.538 

0.536 

0.534 

0.532 

60 

r 

0.468 

0.468 

0.466 

0.465 

0.464 

0.464 

0.464 

0.463 

65 

u.„ 

0.385 

0.385 

0,386 

0.387 

0.389 

0.391 

0.393 

0.395 

70 

0.295 

0.298 

0.299 

0.303 

0.310 

0.315 

0.321 

0.325 

0.329 

75 

0.203 

0.209 

0.214 

0.244 

0.236 

0.247 

0.255 

0.262 

0.268 

80 

0.113 

0.126 

0.136 

0.153 

0.172 

0.186 

0.197 

0.206 

0.213 

85 

0.0361 

0.0610 

0.0751 

0.0969 

0.119 

0.135 

0.148 

0.15* 

0.165 

90 

0 

0.0265 

0.0390 

0.0594 

0.0809 

0.0961 

0.108 

0.117 

0.124 

'awl 


**-/ 


»/2 


RW«)  sin  <p ,.  dip, 


s > 


for  0,  * signal  zenith  angle, 


Also, 


for 


R(0S)  = Sea  Surface  reflectance  as  a function  of  zenith  angle. 
[fl(»s)-f2(0s)]f3(0s) 

s [w  + W]  [W  + W]  ’ 

fj(0s)  = sin2  (0S)  n2  + cos  20s  j , 

f2^s^  = sin  sin  2^s  [ n2  " sin2|^s  ] n ’ 
f3(0s)  = n2  cos20s  - sin20s  cos  20s,  and 
n = sea-water  index  of  refraction. 


(S-8) 


(S-9a) 


(S-9b) 

(S-9c) 

(S-9d) 
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for  v ■ surface  wind  speed  in  knots,  and  the  index  of  refraction  of  sea  water  has  been  taken 
* 4/3.  Equation  (S-l  la)  is  evaluated  in  table  7 and  figure  9 for  v in  knots  (and  m/sec)  and 
M^w  in  ra(ii&ns  (and  degrees). 

The  relative  contribution  of  M\w  to  distribution  of  radiance  at  the  receiver 
will  be  discussed  in  secton  2.3.8.  Except  for  the  clearest  water,  this  has  a small  effect;  and 
so  the  impact  of  neglecting  zenith  angle  effects  and  dissimilar  wave  slopes  in  the  downwind 
and  crosswind  direction  may  be  negligible.  We  therefore  adopt  this  model,  and 

^AW  = °.ropt>10-  (S*1,b) 

until  better  information  is  available. 


Table  7.  Half-angle  rms  air-water  interface  effects 


v.  Wind  Speed 

,W 

Knots 

Metres/Sec 

Milliradians 

Degrees 

0 

0 

0 

0 

2 

1.03 

14.6 

0.84 

4 

2.06 

20.7 

1.18 

8 

4.12 

29.2 

1.67 

14 

7.21 

38.6 

2.21 

20 

10.3 

46.2 

2.65 

26 

13.4 

52.6 

3.0 

32 

16.5 

58.4 

3.35 

38 

19.6 

63.6 

3.64 

SURFACE  WIND  SPEED,  rmtrn/w 


8URFACE  WIND  SPEED  (KNOTS),  V 


Figure  9.  Half-angle  rms  air-water  interface  effects  as  a function  of  wind  speed  v. 
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2.3.6  Relative  Surface  Foam  Coverage 

The  energy  transmission  of  the  air-water  interface  is  composed  of  the  factors: 

raw  =(rawl)X(Taw2),with 

Taw  = total  energy  transmission  of  the  sir-water  interface,  where 

rawl  = a'r*water  interface  transmission  due  to  index  of  refraction  discontinuity 
and 

• aw2  = air-water  interface  transmission  due  tc  foam  and  streaks  on  the  water  surface. 

This  section  treats  raw2,  while  rawl  was  discussed  in  section  2.3.4.  Independent  of 
the  cloud  conditions,  a model  relating  the  fraction  of  surface  covered  to  the  surface  wind 
speed  is  given  by  (ref  8): 

Cf  =<1.2  X 1CT3)  v3-3.  v <9m/sec,  and  (S-12) 


Cf=(1.2(lCT5)  v 3 3 (0.225  v -0.99),  v >9m/sec;  (S-13) 

where, 

Cf  = fraction  of  surface  covered,  and 
v = surface  wind  speed  in  metres/sec. 

Assuming  the  foam  and  streaks  have  an  albedo  of  1 , the  value  of  raw2  is  given  by: 

Taw2  = 1 -0  2 (10*3))  v3'3,  v <9m/sec,  and  (S-!4a) 

Taw2  ® 1 - (!-2  (10_5»  v3,3  (0.225v  - 0.99);  v > 9m/sec  (S-I4b) 

Equations  (S-14a,  b)  are  evaluated  in  table  8 and  figure  10  for  V in  metres/sec  (and  knots). 
Although  this  model  neglects  zenith  angle  effects,  it  is  adopted  pending  further  experimental 
work. 


2.3.7  Water  Energy  Transmission 


The  energy  transmission  of  the  water  portion  of  the  propagation  path  is  denoted  by 
rw.  For  most  water  types  and  receiver  depths,  the  dominant  cause  of  attenuation  is  the 
diffuse  attenuation  coefficient  (ref  8)  of  the  water,  k.  We  therefore  use  the  model: 


rw  = exp  [-  (k  X (path  length  in  water)) . 


(S-15a) 


The  path  length  ir  the  water  is  given  by: 


D/cos  (4>™)  , 


8.  HR  Gordon  and  MM  Jacobi,  Albedo  of  the  Ocean-Atmoipheric  System:  Influence  of  the  Sea  Foam, 
Appl  Optici,  Vol  1^5  (8),  p 2257-2260,  August  1977, 
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where 


D = receiver  depth,  and 
<j>^  = in-water  signal  zenith  angle. 
From  Snell’s  law, 


n sin  = sin  <j>$  , 


where 


n * water  index  of  refraction,  and 
<f>s  “ in-air  signal  zenith  angle. 


Table  9 and  figure  1 1 show  the  values  of  equation  (S-15)  for  n = 1.33.  Since  many  areas  of 
the  ocean  have  a layered  structure  for  k,  we  modify  our  basic  equation  to  yield 


where  we  have  assumed  there  are  j layers  which  differ  in  their  k values,  but  not  in  their 
indices  of  refraction.  Therefore,  the  adopted  model  is: 


(S-l  5b) 


This  model  is  uncertain  in  many  ways: 

(1)  the  value  of  kj  to  use, 

(2)  the  values  of  Dj,  and 

(3)  its  applicability  in  very  clear  water  and/or  at  shallow  receiver  depths. 

Consequently,  although  it  is  the  best  model  available,  it  should  be  revised  when 
better  information  becomes  available. 

2.3.8  Water  Distribution  of  Radiance 

There  is  no  experimentally  verified  expression  for  the  in-water  distribution  of  signal 
radiance,  as  a function  of  incident  beam  collimation,  incident  beam  zenith  angle,  water 
properties,  and  receiver  depth.  Any  model  expression  must  take  into  account  the  following 
items: 


Table  9.  Relation  between  in-air  and  in-water  signal  zenith  angles 
(assuming  sea-water  index  of  refraction,  n * 1 .33) 


<j>t,  In-Air  Signal  Zenith  Angle 
(degrees) 

In-Water  Signal  Zenith  Angle 
(degrees) 

0 

0 

5 

3.75 

10 

7.48 

IS 

11.19 

20 

14.86 

25 

18.48 

30 

22.02 

35 

25.48 

40 

28.82 

45 

32.03 

50 

35.07 

55 

37.91 

60 

40.51 

65 

42.82 

70 

44.81 

75 

46  42 

80 

47.61 

85 

48.34 

Figure  1 1 . Relation  between  in-air  and  in-water  signal  zenith  angles  (assuming  sea-water 
index  of  refraction,  n ■ 1.33). 


( 1 ) the  in-air  zenith  angle  of  the  signal, 

(2)  the  radiance  distribution  incident  on  the  water, 

(3)  the  air-water  interface  effects, 

(4)  the  in-water  zenith  angle  of  the  signal, 

(5)  the  in-water  scattering  effects,  and 

(6)  the  in-water  absorption  effects. 

In  considering  these  last  two,  we  note  that  the  energy  should  decrease  away  from 
zenith  due  to  absorption,  depending  on  depth  and  the  water  properties.  We  have 
therefore  considered  the  following  zenith  angle  dependencies  for  the  radiance  at  the  sub- 
marine receiver: 


(1) 

(2) 


Uniform 


(3) 

(4) 

(5) 


Moreover, 

0W  = angle  measured  from  the  axis,  or  principal  ray  direction,  of  the  in-water  signal 
radiance; 


and 

(J>q  = angle  at  which  the  signal  radiance  reaches  a benchmark. 


We  adopt  (4),  1 - 


sin  0 
sin  0 o 


2 


because  it  is  easy  to  work  with  and  appears  to 


give  a reasonable  representation  of  the  assumed  radiance. 

<Pq  is  related  to  the  half-pov.er  point  of  the  received  radiance  by  the  equation: 


1 - COS  (0i/2> 

1 3 sin^  4>q 

cos  <t>\/2  sin ^ 0i/2  + 2 cos  0y2  - 2 

2 1 - cos  <fiQ  - L 

3 sin^  00 

cos  00  sin'’’  00+2  cos  00  -2 

(S-16a) 
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Equation  (S-l  6a)  is  evaluated  in  table  10.  Values  between  those  shown  are  obtained 
by  linear  interpolation. 


Table  10.  Relation  of  radiance  zero  point,  $q, 
and  received  radiance  half-power  point,  for 
1 - (sin  $w/sin  radiance  distribution 


d>V4  (degrees) 

4>q  (degrees) 

3.8 

5 

7.6 

10 

11.4 

15 

15.2 

20 

19.0 

25 

22.7 

30 

26.5 

35 

30.2 

40 

33.9 

45 

37.5 

50 

41.1 

55 

44.6 

60 

48.1 

65 

51.6 

70 

54.9 

75 

58.2 

80 

We  adopt  the  NOSC  developed  expression  for  the  half-power  point  in  terms  of 
incident  beam  divergence,  air-water  angular  effects  and  in-water  scattering  effects,  and  add 
them  up  as  if  they  were  three  statistically*  independent  effects: 


(S-l  6b) 


for 


fw  = water  contribution 


W 

COS  0S 


> T opt  ' 


(S-l 7a) 


"This  is  actually  an  empirical  result,  and  appears  to  fit  the  NOSC  experimental  results.  A completely 
consistent  theory  of  all  these  effects  would  not  use  the  statistical  independence  as  the  justification 
for  this  expression. 
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Here, 

0£j  = mean  square  single  scattering  angle  in  water, 

s = scattering  coefficient  in  water, 

D = receiver  depth, 

<t> * - in-water  signal  zenith  angle, 


(S-17b) 
(S-l  7c) 

(S-17d) 
(S-l  7e) 


1 
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The  radiance  distribution  enters  into  the  expression  for  the  received  energy  by  the 
expression  r , . -o 


nR  . / sin 


fT® R,  0Q>  ' 


\ L \sin*o/ 


where, 

= solid  angle  of  the  receiver, 

J2q  * full  solid  angle  which  all  the  energy  is  within, 

0^  = half-angle  of  the  receiver  field  of  view,  and 

6 - off-set  angle  between  receiver  optical  axis  and  axis  of  the  incoming  light. 

Therefore,  for  the  general  case, 


f(0R>  <*>0’  6)  * 


$w  = cos'*  | cos  0W  cos  6 + sin  sin  6 sin  0 } 

For  perfect  alignment  between  the  received  light  axis  and  the  receiver  (6  •*  0),  the 
integral  can  be  analytically  evaluated,  with  the  result 

I - cos  !- [cos  0 ^ sin2  0^  + 2 cos  - 2] 

3sin20 


1 - cos  6 q [cos  <Pq  sin-4  + 2 cos  0q  - 2] 

3 sin2  0q  (S-18b) 

Equation  (S-20b)  has  been  evaluated  in  figures  12,  13  for  Qq  = 30  degrees  (figure  1 2) 
and  = ?0  degrees  (figure  13)  and  6 = 0,  20,  30,  40,  50,  and  60  degrees. 

These  results  will  be  used  in  section  2.3. 10. 


2.3.9  Received  Pulse  Width/Shape 

There  is  no  verified  model  which  predicts  the  received  pulse  width  and  shape  at  the 
submarine  receiver  as  a function  of  cloud  properties,  cloud  height  above  the  water,  water 
properties,  and  the  receiver  field-of-view.  It  is  not  intended  to  develop  such  a final  model 
now,  but  only  to  present  some  reasonable  expressions  needed  for  the  provisional  model  here. 

It  is  the  universal  experience  of  experimenters  (both  empirical  and  computer  simula- 
tion) that  after  penetration  of  a multiple  scattering  medium,  the  received  pulse  has  a short 
rise  time  and  a long  falltime.  Because  of  its  nearness  to  Bucher's  computer  simulations 
(ref.  2,  as  shown  in  figure  14)  we  shall  assume  the  underwater  receiver  sees  a pulse  shape 
given  by: 

f(t)  = te“kt  . (S-19a) 


The  properties  of  this  form  are  as  follows: 

maximum  value  of  f(t)  occurs  at  tm  = ~ , (S-19b) 

maximum  value  of  f(t)  = tm  e~*  * 0.368  tm  (S-19c) 

average  value  of  t,  defined  as 


> 


= 2 t 


m • 


(S-20a) 


Half-power  points  of  f(t)  occur  at  0.233  tm  and  2.68  so  the  width  between 
half-power  points  = 2.45  tm.  * ’ (S-20b) 

The  rms  value  of  t,  defined  as 


Hence, 


and  its 


f t2  fit)  d t 

“A 


CO 

f f(t) 

J A 


d t 


Vi 


= s/6  t 


m 


the  variance  of  t is  given  by 
var  t = (t2  -T2)  = 2 t2, 
standard  deviation,  by 
-o  yi 


rt2-t2r=y2  * 


m 


(S-20c) 


(S-21a) 


(S-21b) 


Figure  14.  Comparison  of  f(t)  and  Bucher’s  Monte  Carlo  pulse  shape. 

The  area  under  the  f(t)  curve  is  given  by 

OO 

[ f(t)dt  = t^  . (S-22) 

Jo 

In  principle,  there  are  three  significant  contributions  to  the  received  pulse  width, 
which  we  shall  define  as  the  time  between  half-power  points, 

(S-23) 

where  we  have  neglected  the  initial  pulse  width  at  the  transmitter,  assumed  that  the  effects 
add  serially,  not  statistically,  and  where 

Atw  » pulse  width  due  to  water  portion  of  the  path, 

Atcw  * additional  pulse  width  due  to  cloud-to-water  path,  and 

Atc  = pulse  width  after  emerging  from  the  cloud. 

The  additional  pulse  width  due  to  the  water  portion  of  the  path  is  caused  by  the 
in-water  multiple  scattering,  and  so  it  occurs  in  the  absence  or  presence  of  clouds.  As  seen 
in  figure  1 5,  a signal  at  a zenith  angle  of  0 degrees  and  a receiver  of  half-angle  d ^ at  a 
depth  D leads  to  a pulse  width  (for  a uniform  contribution  throughout  that  field  of  view),  of 


(S-24) 


i 

i 

f 

l 
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For  thin  clouds  (ropt  < 10)  we  have  found  no  verified  expression  for  the  cloud 
effects  or  the  cloud  to  water  effects.  Since  equations  (S-28)  and  (S-29)  already  imply 
stretching  to  a few  hundred  nanoseconds,  we  can  neglect  the  thin  cloud  effects,  and  take 

(S-26a) 


For  thick  clouds  we  adopt  the  Stotts  (ref.  9)  expression,  so  that: 


(S-26b) 


where, 

c = speed  of  light  in  air 

T = cloud  geometrical  thickness, 

cjq  = single  scatter  albedo,  and 

6 = mean  scattering  angle  in  the  cloud. 

Equation  (S-26b)  has  been  evaluated  in  table  1 I and  figure  17  for  the  typical  values 
of  w0  * 1 , d = 0.66  rad  (37  degrees),  and  fixing  T * ropt/oc  * 25  r0pt  suitable  for  a strato- 
cumulus  cloud.  (Also  shown  for  comparison,  is  the  result  estimated  in  ref  2.) 

These  values  probably  overestimate  the  pulse  widths  at  the  lower  values  of  ropt,  but 
they  are  adopted  until  a verified  model  for  all  vaues  of  ropt  is  developed. 

We  use  the  values  from  table  1 1 for  ropt  = 20,  40,  and  60  along  with  the  normalized 
pulse  shape  (f(t)  = (tj^)  t exp  - (t/tm)l , to  plot  representative  pulses  in  figure  18.  The 
drastic  dependence  of  the  pulse  height  and  width  on  optical  thickness  for  the  assumed  mode) 
is  clearly  seen. 

A careful  analysis  of  the  thick  cloud  to  water  propagation  region  has  found  that  it 
adds  a negligible  amount  relative  to  Atc  for  reasonable  cloud  properties  and  cloud  heights. 
Hence,  we  take 

(S-25b) 


NORMALIZED  PULSE 


T 1 | 1 — I 1 1 1 I 1 I I T 


Figure  J8.  Representative  normalized  pulse  shapes  as  a function  of  cloud  optical  thickness. 
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2.3.10  Overall  Signal  Equations 


The  optical  detection  mechanism  responds  to  the  received  energy  as  a function  of 
time,  i.e.,  the  instantaneous  optical  power.  The  total  received  optical  energy  and  the 
instantaneous  optical  power  are  related  by: 


(S-27) 


ER  = total  received  optical  energy  per  pulse,  and 


PR(t)  » instantaneous  received  optical  power. 

Writing 

PR(t)  = Ae  f(t), 

(S-28) 

for 

AE  = normalization  parameter,  and 

f(t)  = received  optical  pulse  shape, 

then 

. _ er 

AE-  oo 

(S-29) 

f f(Odt 

Jo 

so  that 

Er  f(t) 

(S-30) 

pr(0=  — . 

f f(t)  dt 

Jo 

Using  the  received  pulse  shape  of  section  2,3.9, 

f(t)  = t exp  - p~  , 

(S-3 1 ) 

(S-32) 

and 


pR(t)  = Er 


exp  - (t/tm) 


(S  -33) 
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for 


tm  = (2.45)“*  X (time  between  half  power  points). 


(S-34) 


The  total  received  optical  energy  is  given  by  the  range  equation: 

Er  * (transmitted  energy  per  pulse)  X (transmitter  optics  transmission) 

^ / solid  angle  subtended  by  the  receiver  at  the  source  j 
V solid  angle  transmitted  into  > 

X (clear  atmospheric  energy  transmission) 

X (cloud  energy  transmission) 

X (cloud  to  water  energy  transmission ) 

X (air-water  interface  energy  transmission) 

X (water  energy  transmission ) 

X ( receiver  optics  transmission ) 

X (fraction  of  incident  radiance  within  receiver  field-of-view  ) • 


(S-35) 


We  take 


tt 

tr 

d 

R 

<*d2/4)  _ 


>2 


6~ 


transmitted  energy  per  pulse, 
transmitter  optics  transmission, 
receiver  optics  transmission, 

diameter  of  receiver  aperture, 
range  from  source  to  receiver, 

solid  angle  subtended  by  the  receiver  at  the  source, 

full  angle  e-2  irradiance  angle  transmitted  into*, 


ir6  2T 


solid  angle  transmitted  into, 

= clear  atmosphere  energy  transmission  (as  discussed  in  section  2.3.1), 

= cloud  energy  transmission  (as  discussed  in  section  2.3.2), 

= cloud  to  water  energy  transmission  (as  discussed  in  section  2.3.3), 

= air-water  interface  energy  transmission  (as  discussed  in  sections  2.3.4 
and  2.3.6), 

* water  energy  transmission  (as  discussed  in  section  2.3.7),  and 
I(0W)  = water  radiance  distribution  (as  discussed  in  section  2.3.8). 

Therefore,  the  received  optical  energy  is  given  by 


' cw 


'aw 


' w 


Ed  = 


^Ep7T 


r(*d2/4)  yR\ 

) 


ra  Tc  rcw  Taw  rw  ^o*  ^R* 


(S-36) 


•This  assumes  such  large  spots  that  additional  cloud  and  water  spreading  is  negligible. 
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The  fraction  of  the  incident  radiance  within  the  receiver  field-of-view  is  given  by  (for 
perfect  alignment  between  beam  axis  and  receiver  axis): 


I ($w)  d Si 

. ; 

I ($w)  d SI 

for 

* half  angle  of  the  receiver  field-of-v.iew,  and 
<t>0  a off-axis  angle  at  which  incoming  radiance  equals  zero. 
Using  the  model  adopted  in  section  2.3.8, 


(S-37) 


f($2)  * f(#Q>  0r) 


1 - cos  - 1- — I cos  sin2  0^  + 2 cos  0^  - 2 I 

3sin2$0  1 1 

cos  $q  sin2  $o  + 2 cos  - 2 J 


(S-38a) 


1 - cos  ( 


9 

3 sin  $o 


and,  f($o,0£)=l,  for0£>$o  . (S-38b) 

Using  equations  (S-36),  (S-27),  and  (S-33)  results  in  the  evaluation  of  the  instantaneous 
optical  power  in  terms  of  all  the  other  signal  group  sub-models. 
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3.0  NOISE  PROPAGATION  GROUP 


This  section  discusses  the  proposed  equation  group  for  the  propagation  of  the  noise 
relative  to  a single  signal  pulse. 

3.1  OBJECTIVE,  STRUCTURE,  AND  FUNCTION 


s~  | This  section  considers  the  basic  approach  used  in  the  detailed  noise  calculation 

k i modules  presented  in  section  3.3,  and  presents  flow  charts  showing  the  interrelationship  of 

- j ' these  sub-models  and  their  required  inputs.  (These  inputs  are  discussed  in  more  detail  in 

:■ 1 ; section  3.2.) 

i 

; | i 

| | 3.1.1  Philosophy  of  Approach 

Since  all  the  background  sources  that  must  propagate  through  clouds  are  continuous 
' in  time,  that  part  of  the  modeling  related  to  the  temporal  effects  is  not  present  here.  This 

f simplifies  the  noise  modeling. 


' i 
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i 


ti 
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On  the  other  hand,  the  sources  of  the  noise  are  so  different  (sun,  moon,  sky-light, 
star-light,  bioluminescence,  shot  noise  in  the  receiver  amplifier,  detector  dark  current,  etc.), 
that  is  has  been  difficult  to  develop  a single  unified  approach  to  all  of  them.  The  approach, 
therefore,  is  described  below. 

Express  the  noise  contribution  as  a noise-equivalent-optical-power  (NEPtot),  which 
is  the  root  sum  of  squares  of  all  the  noise-equivalent-optical  power  contributions  of  the 
individual  noise  components.  The  NEP{0t  is  then  directly  comparable  to  the  instantaneous 
received  optical  power  Pj^(t)  developed  in  section  2,  and  the  signal-to-noise  ratio  is  PR/NEPtot; 
for  Pr  the  peak  value  of  the  received  signal  power. 

Take  the  out-of-water  background  sources  in  terms  of  equivalent  exo-atmospheric 
radiances,  and  then  their  propagation  through  the  atmosphere  and  water  path  is  treated  in 
parallel  with  the  signal  energy  transmission  of  section  2.3-  The  angular  effects  and  noise 
radiance  distribution  are  also  treated  in  a manner  similar  to  that  of  the  signal  in  section  2. 

Take  the  noise  contributions  of  the  background  sources  as  the  1 -sigma  point  in  the 
fluctuations  generated  in  the  signal  current  by  their  steady  presence,  and  express  their 
contributions  in  terms  of  an  equivalent  optical  power. 

Treat  the  amplifier  shot  noise,  detector  dark  current  and  signal  shot  noise  in  the 
standard  way,  and  express  their  contributions  in  terms  of  an  equivalent  optical  power. 


Present  the  modeling  in  modular  fashion,  so  that  the  effect  of  each  portion  of  the 
path  is  evident.  In  addition,  as  further  experiments  and  analysis  are  undertaken,  pieces  of 
the  model  may  be  upgraded  without  requiring  extensive  modifications  to  this  publication 
or  the  total  model. 


Separate  the  cloud  conditions  into  clear/thin  cloud,  corresponding  to  an  optical 
i thickness  (r opt  ) < 1 0 and  thick  cloud  for  r opt  >10.  Below  r 0 pt  = 1 0 one  set  of 

noise-modules  is  assumed  to  apply,  while  above  it  a different  set  applies.  In  many  cases 
1 these  noise  modules  do  not  correspond  at  r0p(  = 10,  and  so  the  overall  model  should 

only  be  used  for  rQpt  <10androp{  >10.  (Further  analysis  and  experiments  on  the 
“multiple  forward  scattering”  region  should  enable  these  sub-models  to  be  upgraded, 
and  this  inconsistency  removed.) 
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Finally,  we  assume  appropriate  and  simple  analytic  forms  for  at-present  unknown 
functions  such  as  the  radiance  distributions.  This  enables  us  to  generate  analytic  results 
(except  for  the  receiver  axis  offset  from  the  beam  axis  of  the  incident  radiation),  which 
is  an  aid  to  a physical  understanding  of  the  situation. 

3.1.2  Noise  Group  Structure 

A schematic  of  the  overall  noise  equivalent  power  propagation  group  is  shown  in 
figure  19.  Given  the  input  parameters,  the  path  energy  transmission  and  angular  and 
radiance  distributions  are  derived  for  the  four  “exo-atmospheric”  background  sources. 
Then,  using  additional  input  data,  the  average  background  power  for  all  background  sources 
is  derived,  and  so,  the  noise -equivalent  optical  power  for  all  the  noise  components. 


Figure  19.  Schematic  of  single-pulse  noise  equivalent  power  downlink  propagation  group. 

Figure  20  is  a detailed  flow  diagram  of  the  direct  sunlight  contribution,  showing  the 
calculations  that  must  be  performed  to  arrive  at  the  required  output.  (The  flow  charts  for 
the  moonlight,  blue  sky-light  and  starlight/zodiacal  light  are  identical  to  this  one,  while  the 
simpler  one  for  the  bioluminscence  is  shown  in  figure  2 1 .)  The  following  six  steps  are  used 
to  exercise  the  noise  propagation  group: 

• The  input  parameters  are  listed  in  the  eight  ellipses  on  the  left-hand  side  of  the 
figure  including  source,  clear  atmosphere,  cloud,  cloud  to  water,  air-water  inter- 
face, water,  receiver  and  signal  characteristic  parameters.  (The  symbols  are  defined 
in  the  glossary  and  also  in  the  input  discussion  in  section  3.2.) 

• The  primary  calculation  equations  are  represented  by  the  rectangular  boxes. 

Within  each  box  is  the  symbol  for  the  parameter  to  be  calculated  and  the 
equation  number  (from  section  3.3)  for  the  equation  to  be  used  to  calculate 
that  parameter. 

• The  first  quantity  calculated  is  the  cloud  optical  thickness,  r0p(,  since  this  deter- 
mines which  equations  should  be  used  to  calculate  many  other  parameters.  Whenever 
//  appears  in  a rectangular  box,  the  equation  number  preceeding  it  refers  to 

t opt  > 10.  while  the  equation  number  following  it  refers  to  r 0pt  < 10.  Hence, 
given  the  value  of  r opt,  the  rest  of  the  models  to  be  used  are  specifically  determined. 
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LEGEND: 


PARAMETER: 

THICK  CLOUD  E'Q'N// 
THIN  CLOUD  E'Q'N 


PARAMETERS  DEFINED  IN  GLOSSARY 
INPUT  DISCUSSED  IN  SECTION  3.2 
EQUATIONS  DEVELOPED  IN  SECTION  3.3 


THE  FLOW  DIAGRAM  FOR  P_„,  Pfte,  AND  P, 
mu  do  c 

IS  EQUIVALENT  TO  THIS  ONE  THROUGH  P^ 
Pe(_  FLOW  DIAGRAM  IS  SHOWN  IN  FIGURE  4-C 


Figure  20.  Flow  diagram  of  tingle-pulae  downlink  noiae  propagation  group. 
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• The  second  set  of  calculations  performed  are  of  two  types: 

(a)  path  transmission,  including  r'a,  t'c,  t'cw,  raw,  and  r^,;  and, 

(b)  angular  and  radiance  distribution  including  fa,  faw,  fw,  0O  and  f'  (0g  4>0,  6). 

• The  path  transmission,  angular  and  radiance  distribution,  source  and  receiver  para- 
meters are  then  used  to  calculate  the  average  background  power  due  to  that  source. 

• The  total  average  background  power  due  to  all  sources  is  then  calculated. 

• The  total  average  background  power,  receiver  and  signal  characteristics  are  then  used 
to  calculate  the  total  noise  equivalent  optical  power. 

Figure  21  shows  the  flow  chart  for  calculating  the  average  background  power  due  to 
bioluminescence,  Pgj^.  This  value  of  Pgj_  enters  into  figure  20  in  the  2g  calculations. 


SOURCE 


RECEIVER 

tfR;7R;d; 


f>0L  : (N-S1) 


TO  2 Pg  (N-6?) 


LEGEND 


PARAMETER: 

EQUATION 


PARAMETERS  DEFINED  IN  GLOSSARY 
INPUT  OISCUSSlON  IN  SECTION  3.2 
EQUATIONS  DEVELOPED  IN  SECTION  3.3 


Figure  21 . Flow  diagram  of  single-pulse  downlink  noise  propagation 
submodel  (noise  due  to  bioluminescence). 

3.2  INPUT  FOR  NOISE  CALCULATION 

This  section  discusses  the  form  (and  values  in  some  cases)  of  the  required  inputs  to 
the  noise  group,  in  terms  of  the  eight  categories:  source,  clear  atmosphere,  cloud,  cloud  to 
water,  air-water  interface,  water,  receiver  and  signal  characteristics. 

3.2.1.  Source 

There  are  five  sources  of  the  average  background,  which  are  treated  here  as  inde- 
pendent. These  sources  are:  sunlight,  moonlight,  blue  skylight,  starlight/zodiacal  light, 
and  bioluminescence. 

It  is  expected  that  these  sources  will  be  treated  rationally  when  using  them  as  inputs, 
so  that  when  sunlight  is  present,  only  the  blue-skylight  will  be  expected  to  be  a non-zero 
value;  and  when  moonlight  is  present  only  the  starlight/zodiacal  light  and  bioluminescence 
will  be  used  with  non-zero  values.  (A  full  diurnal  model  for  these  sources  is  not  considered 
appropriate  at  this  time.) 
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We  consider  each  of  the  five  noise  sources  separately. 
SUN 


Symbol 


Description 


Units 


6%n 


H 


sX 


The  half-angle  subtended  by  the  sun  at  the  earth.  Its 
value  is  taken  as  (ref  10):  4.65  10“^  radians. 


Effective  exo-atmospheric  spectral  radiance  of  the  sun. 
The  value  assumed  will  be  2000  watts/m^  - p. 


radians 

watts/( metres)^ 
(microns) 


MOON 


Symbol 


Description 


Units 


*m/2 

^muX 


The  half-angle  subtended  by  the  moon  at  the  earth.  radians 

Its  value  is  equal  to  (ref  1 1):  6^2  = 4.65  (10“^)  radians. 

watts 

Effective  exo-atmospheric  spectral  irradiance  of  the  r 

moon.  It  will  be  assumed  to  equal  4.3  watts/m^  - p.  (metres)  (microns) 


STARLIGHT/ZODIACAL  LIGHT 

Symbol  Description 

Lz  Effective  exo-atmospheric  spectral  radiance  due  to  all 
non-lunar  night  time  sources.  The  value  (ref  1 2)  is 
3(10~6)  (watts/(metres)2(sradXmicron)l  in  the  blue- 
green  spectral  region. 


Units 


(watts/(metres)^ 
(steradians) 
(microns)  1 


BIOLUMINESCENCE 

Symbol  Description 

LjjL  Spectral  irradiance  of  the  ambient  bioluminescence 

sources  at  the  aperture  of  the  submarine  receiver.  The 
value  for  this  parameter  is  least  well  known  of  all  the 
background  contributors.  As  previously  done,  and  until 
further  information  is  available,  we  use  the  values 
provided  in  the  SAOCS  RFP,  so  that  (ref  1 3)  Lg^ 

~ ( 10“^)  watts/m^  microns. 


Units 

(watts/(  metres)^ 
(microns)) 


BLUE  SKYLIGHT 

Symbol  Description  Units 

lb  Effective  exo-atmospheric  spectral  radiance  of  [watts/fmetres)^ 

the  blue-sky.  Its  value  is  assumed  to  be  250  (microns)) 

watts/{(metres)2(srad)(micron>! 


10.  Handbook  of  Geophysics  (rev  ed),  p 1 7-1 , 7-2.  MacMillan  Co,  New  York,  1960. 

1 1.  RC  Haynes,  Introduction  to  Space  Service,  p 4-5,  John  Wiley  and  Sons,  New  York,  1971. 

12.  WK  Pratt,  Laser  Communication  Systems,  ref  4,  p 122,  fig  6-7,  John  Wiley  and  Sons,  New  York,  1969. 

13.  NOSC  Interim  Report  No.  2,  Contract  No  N00039-77-C-0100,  Optical  Submarine  Communication  by 
Aerospace  Relay  (OSCAR)  (U),  T Flom.PFTitterton.et  al.,p3-22to  3-24,  Secret,  1 May  1978. 
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3.2.2  Clear  Atmosphere 

i 

i 

; 

The  required  parameters  are: 

4 

c 

* 

r 

Symbol 

Description 

Units  ^ 

- 

b 

Effective  clear  atmospheric  optical  thickness.  For  a 

■i 

i 

r 

zenith  transmission  of  70%;  b = 0.357 

1 

i 

^su 

In-air  solar  zenith  angle 

radians  ] 

l 

^mu 

In-air  lunar  zenith  angle 

a 

radians 

t 

1 

.4 

\ 

3.2.3  Cloud 

a 

i 

if.  • 

The  required  parameters  are: 

1 

j 

] 

Units  | 

• 

Symbol 

Description 

; 

T 

Geometric  or  physical  thickness  of  the  cloud. 

metres  a 

°c 

Average  extinction  coefficient  of  the  cloud. 

1 1 

(metres)-  ? 

. i 

^su 

In-air  solar  zenith  angle 

4 

radians  1 

^mu 

In-air  lunar  zenith  angle 

radians  1 

* 

-f 

<cos  d> 

The  average  value  of  the  cosine  of  the  scattering  angle 

i 

for  single  scattering  within  the  cloud. 

tj 

3 

3 

3.2.4  Cloud  to  Water 

1 

■i 

i 

No  parameters  in  this  area  affect  the  noise  properties. 

j 

* 

3.2.5  Water 

j 

The  required  parameters  are: 

1 

i 

Symbol 

Description 

Units  1 

ki 

Diffuse  attenuation  coefficient  of  the  i’th  water  layer. 

(metres)-'  j 

Di 

Thickness  of  the  i’th  water  layer 

(metres)  1 

c 

In-water  solar  zenith  angle 

i 

radians  j 

3 

• 

.w 

^mu 

In-water  lunar  zenith  angle 

radians  | 

n 

Water  index  of  refraction 

• 

^si 

Root-Mean-Square  angle  for  a single  scattering  event  in 
the  water. 

radians  * 

s 

Scattering  coefficient  for  the  entire  water  path 

(metres)-'  | 

D 

Depth  of  the  submarine  receiver 

metres  j 

i 

• 

S 

fc. 

V 

- *=* 
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3.2.6  Air/Water  Interface 


The  required  parameters  are: 


Symbol  Description 

0SU  In-air  solar  zenith  angle 

0mu  In-air  lunar  zenith  angle 

v Surface  wind  speed 

n Water  index  of  refraction 

3.2.7  Receiver 


The  required  parameters  are: 


Symbol 

6 

TR 

d 

Bopt 

(kT) 

Fa 

G 

(r?e/hu) 

RL 

Id 

F 

D 


Description 

Half-angle  of  the  receiver  field  of  view 

Off-set  angle  between  the  in-water  noise  source  beam 
and  the  receiver  optical  axis 

Transmission  of  the  receiver  optical  aperture 

Diameter  of  the  receiver  optical  aperture 
Passband  of  the  optical  filter 

Thermal  noise  contribution  in  the  amplifier 
Excess  amplifier  noise  o . thermal  noise 

Gain  of  the  photo-detector 

Responsivity  of  the  photo-surface 

Load  resistance  following  the  photo-detector 

Dark  current  at  the  detector  cathode 

Excess  noise  in  the  photo-detector  gain 
Depth  of  the  submarine  receiver 


Units 

radians 

radians 

metres/second 


Units 

radians 

radians 


metres 

microns 

Joules 


amps/Wat  is 
Ohms 
amps 

metres 


3.2.8  Signal  Characteristics 

Two  parameters  from  the  signal  characteristics  which  enter  into  the  Noise 
calculations  are: 

Symbol  Description  Unit 

tm  Time  after  pulse  initiation  at  which  it  peaks,  for  a seconds 

t exp  -(t/tm)  shape. 

f^R  Peak  value  of  the  received  optical  signal  power.  watts 
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3.3  NOISE  CALCULATION  MODULES 


This  section  develops  all  the  equations  used  in  the  calculation  of  the  noise  contri- 
bution to  the  total  noise  equivalent  optical  power. 

Sections  3.3. 1 , 3.3.2,  3.3.3,  3.3.4,  3.3.6,  and  3.3.7  consider  the  path  transmission  of 
the  energy. 

Sections  3.3.5  and  3.3.8  consider  the  angular  effects  and  the  distribution  of  the 
received  radiance. 

Section  3.3.9  considers  the  electrical  detection  bandwidth  in  terms  of  the  received 
pulse  width. 

In  each  of  these  sections,  after  the  equations  are  developed  they  are  evaluated  for 
typical  cases  in  both  tables  and  figures. 

Sections  3.3.10  and  3.3. 14  then  derive  tht  average  background  optical  power  for 
the  five  types  of  background  sources,  and  Section  3.3. 1 5 presents  the  expression  for  the 
total  noise  equivalent  optical  power  due  to  all  noise  sources. 

3.3. 1 Clear  Atmospheric  Transmission 

In  the  absence  of  any  clouds  or  aerosols,  the  clear  atmospheric  transmission  is 
described  by  the  term  r'a.  Using  the  approximate  AFCRL  model  (ref  1 ),  the  solar  (or  lunar) 
zenith  angle  dependence  is  given  by: 


r'a  = exp  (-b  sec  0$u)  , (N-la) 


ra  = exp  (-b  sec  0mu) 


(N-lb) 


r'a  - solar  (or  lunar)  clear  atmospheric  transmission 

b = effective  clear  atmosphere  optical  thickness 
0SU  = solar  zenith  angle.  While  conversely, 

0SU  is  replaced  by  0mu  = lunar  zenith  angle. 

The  typical  value  of  b is  determined  by 
r'a(0su  = O)  = O.7  = cxp(-b)  , 


b = 0.357. 

Table  1 2 and  figure  22  show  the  values  of  ra  as  a function  of  solar  zenith  angle. 

The  other  two  out-of-water  sources  of  background  radiation  are  taken  as  uniformly 
distributed  over  the  hemisphere.  Then  the  effective  atmospheric  transmission  is  weighted 
by  the  transmission  at  each  zenith  angle,  or,  for  the  blue-sky  and  the  stellar  sources, 


£71  <1/  £ 

f d0  f sin  9 d6  exp  - (b  sec  0] 

J o J o 

2ir  ?r/2 

f d0  sin  6 d0 

J o o 


53 


3.3.2  Cloud  Energy  Transmission 

The  transmission  of  sunlight  or  moonlight  by  clouds  does  not  have  an  experimentally 
verified  expression.  Most  work  in  cloud  transmission  has  been  a broadband  treatment  of 
transmitted  irradiance  (watts/m^),  and  there  have  been  no  experiments  which  collected  the 
total  energy  emerging  from  thick  or  thin  clouds. 

We  propose  to  adopt  a number  of  different  (but  consistent)  models,  depending  on 
the  characteristics  of  the  noise  source. 

For  the  sun  and  moon  we  shall  adopt  the  same  model  as  for  the  signal  (section 
2.3.2),  with  the  exception  that  both  the  sunlight  and  moonlight  do  spread  out  with  zenith 
angle,  and  thus  the  extra  cosine  factor  is  always  present. 

Therefore,  for  the  sun, 


(N-2a) 


(N-2b) 


(N-2c) 


where 

r'cs  = cloud  energy  transmission  of  direct  sunlight, 
ropt  = optical  thickness  uf  the  cloud 

<cos  d>  = mean  cosine  of  the  scattering  angle, 
and  0SU  = solar  zenith  angle.  A(^s)  was  given  in  equation  (S-4)  of  section  2.3.2. 
As  before, 


Tcs: 


A(#su) 


ropt(l-<cosf?>)+1.42 


V3O-<cosd»(l-co0) 


r + 142 

°»*  1-  <cos0> 


exp 


\J 3(  1 -<cos0»(l  -uj0) 


1.42 


°Pl  1 - <cos  Q> 


1 - exp  - 


2^3(1  -<cosd»(l  - w0) 


T + _ 1.42 

°Pf  1 - <cos  6> 


X [ cos^su 


1.0 


forTopt>  10- 


Tcs  * 

1 0 085  r f 1 '69  _ 1 1 nc  rh  1 ^ 

0085  ropt  1,0(1  _<cos0»+ 1.421  * 

for^pt4  10- 

rcs  = cos*su-  for ’’opt30 

(N-3) 


Here, 

T - geometrical  thickness  of  the  cloud,  and 
aQ  = mean  extinction  coefficient  of  the  cloud. 
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For  the  moon, 


(N4a) 


(N-4b) 


(N~4c) 

Here, 

r'cm  » cloud  energy  transmission  of  direct  moonlight, 

and 

0mu  = ,unar  zeni,h  angle. 

The  typical  energy  transmission  (for  <cos  6>  = 0.83)  described  by  equations  (N-2) 
and  (N-4)  is  shown  in  table  13  and  Figure  23. 

The  zenith  angle  dependence  for  both  regimes  of  r0pt  is  shown  in  tabic  14  and 
figure  24. 

Again,  there  is  a discontinuity  in  the  zenith  angle  dependence  at  ropt  - 10,  which  we 
shall  ignore  pending  a verified  model  of  cloud  energy  transmission 

The  other  two  out-of-water  sources  of  the  ambient  background  are  approximated 
as  uniform  sources  distributed  across  the  full  hemisphere.  Then,  an  extra  factor  arises  due 
to  their  distribution  in  angle  of  incidence  upon  the  cloud.  This  extra  factor  is  given  by 

r»/2 

/ f(0)  sin  0 d 0 ; 

J o 

for  f(<p)  - A (<f>)  for  thick  clouds, 

= cos  <t>  for  thin  clouds,  and 
<t>  = zenith  angle  of  the  incident  light. 


f1  S 

'em 


I L 

ropt<i"<cosfl>>+,-42i  r 


exp  - 


2^/3(  1 - <cos  0>X  1 - coq) 
1.42 


✓3o-<c«*>xi-«o)  '.pt* 


1 - exp- 


2v/3(  1 - <cos  fl>X  1 - wq)  I f0pt  + 


+ 1.42 

T°Pl  1-<cos0> 


costf 


mu 


1.0 


forropt>  10- 


and 


rcm 


1 - 0.085  r 


opt 


1.69 


Also, 


Tcm  cos  ^mu 


10(1  — <cos  d>)  + 1.42 
forropt<  I0' 


for  ropt  = 0 


(c°s4»miJ>‘  ; 
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NORMALIZED  CLOUD  TRANSMISSION.  T* 


Table  14.  Zenith  angle  dependence  of  sun  and  moon  cloud  energy 
transmission  (normalized  to  dsu  = 0 and  ■ 0). 


‘7  TS  .TJfJR'.Ti  ' 


*su'*mu  Angle 

Thick  Cloud  Dependence 

(v>,0) 

Thin  Cloud  Dependence 

(v<10) 

0 

1 

1 

10 

0.96 

0.97 

20 

0.90 

0.88 

30 

0.79 

0.75 

40 

0.66 

0.59 

50 

0.50 

0.41 

60 

0.34 

0.25 

70 

0.20 

0.12 

B0 

0.08 

0.03 

85 

0.03 

0.008 

Figure  24.  Thin  and  thick  cloud  zenith  angle  dependence  of  cloud  transmission  normalized 
to  zenith. 


Evaluating  this  integral  we  find  it  equals  0.669S  (expressed  2/3)  for  thick  clouds  and 
0.5  (expressed  1/2)  for  thin  ones.  Therefore,  for  the  blue  sky  background,  the  cloud 
energy  transmission  is  given  by: 


1.69 


(1  -<cos0»+  1.42 


1 - <cos  6> 


so -<«.*»( i -«„>  |rop»*  i .~<co2.a> 


rCB  = ",  1 - 0 (185  'opt 

1.69 

10(1  - <cos  0>)+  1.4? 

for  t{ 

t 

>pt<»°. 

rCB  = 

l,forropt  = ° 

The  discontinuity  at  rQpt  = 10  is  again  present,  and  again  neglected  until  a better 
experimental  model  is  derived. 

Note  equation  (N  5a,  b)  are  independent  of  solar  zenith  angle.  However,  the 
strength  of  the  radiance  incident  from  the  blue  sky  does  depend  on  solar  zenith  angle. 

For  the  moonless  night  case,  the  stellar/zodiacal  light  cloud  energy  transmission 
is  given  by 


1.69 


( 1 - <cos  u>)  + 1 .42 


/•n/3(!  - <cos  0>)(  I - u>q) 


1 - <cos  d> 


1 - <cos  9> 


» 


*CZ  = i 1- 0.085  r 


1_69 

°P‘  10(1 -<cosfl»+ 1.42 


forropt<  ,0> 


(N-6b) 


TCZ  = 1*forTopt  = 0 ■ 


(N-6c) 


Table  1 5 and  figure  25  show  a typical  cloud  energy  transmission  as  a function  of 
cloud  optical  thickness  for  these  uniform  sources  of  ambient  background. 

Table  15.  Typical  cloud  energy  transmission  for  blue  skylight 
and  stellar^odiacal  light  (<cos  6>  = 0.83,  cJq  = 1). 

r , , Optical  Thickness  I r' h , r'  Cloud  Energy  Transmission 


027  (thin);  0.36  (thick) 


1 

f 

l 

\ 

I 


* 

1 


r 


*. 


) 


i 
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Figure  25.  Typical  cloud  energy  transmission  for  blue  skylight  and 
stellar  light  (<cos  8>  = 0.83,  Wq  = 1). 

3.3.3  Goud  to  Water  Energy  Transmission 

Because  the  sun  and  rnoon  are  effectively  sources  of  infinite  plane  waves,  and  the 
blue  sky  and  stellar  background  cover  the  entire  hemisphere,  there  is  r.o  “spot”  or  “beam” 
enlargement  in  propagation  from  the  cloud  base  to  the  water  surface.  Therefore,  the 
transmission  of  noise  energy  from  cloud  base  to  the  water  surface  is  given  by: 


(N-7) 


for  all  cloud  conditions. 

3.3.4  Air- Water  Transmission 

The  energy  transmission  of  the  air-water  interface  is  composed  of  two  factors: 
Taw  ~ (rawl ) * (raw2^  ’ 


(N-8) 


for, 


raw  = total  energy  transmission  of  air-water  interface, 

Tawl  = a*r*water  interface  transmission  due  to  index  of  refraction  discontinuity,  and 

Taw2  = a'r*water  interface  transmission  due  to  foam  and  streaks  on  the  sea  surface. 

This  section  treats  r'aw[.  while  r'aw2  is  discussed  in  section  3.3.6.  For  thin  clouds  and  clear 
weather  (r0pt  < 1 0)  the  solar  and  lunar  energy  transmission  is  again  given  as  a function  of 
wind  speed  and  solar/lunar  zenith  angle  in  table  16  and  figure  26. 
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air-sea  interface  transmission, 


Table  16.  T'awls/fawlm  t*me  averted  downlink  air-sea  interface  transmittance 
(for  thin  clouds,  r0pt  < 10). 


Vj  Wind  Speed 


0 

1.03 

2.06 

4 12 

7.21 

10.3 

13.4 

16.5 

0 

2 

4 

8 

14 

20 

26 

32 

0 

0.979 

0.977 

0.976 

0.974 

0.970 

5 

0.975 

0.974 

0,972 

0.970 

0.966 

10 

0.964 

0.962 

0,961 

0.959 

0.955 

15 

0.945 

0.944 

0.943 

0.940 

0.936 

20 

0.920 

0.918 

0.917 

0.914 

0.910 

25 

0.887 

0.885 

0.884 

0.881 

0.877 

30 

0.847 

0.845 

0.844 

0.841 

0.837 

36 

0.800 

0.798 

0.797 

0.794 

0.790 

40 

0.747 

0.745 

0.743 

0.741 

0.736 

45 

0.687 

0.685 

0.684 

0.681 

0.677 

50 

0.620 

0.619 

0.617 

0.615 

0.611 

55 

0.548 

0.546 

0.545 

0.543 

0.540 

60 

0.469 

0.468 

0.468 

0.466 

0.465 

65 
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0.385 

0.385 

0.386 

0.387 
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0.310 

75 

0.203 

0.209 

0.214 
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0.172 

85 

0.0361 
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0.0969 

0.119 

90 

0 

0.0265 

0.0390 

0.0594 

0.0809 

0.967 

0.963 

0.951 

0.933 

0.907 

0.873 

0.833 

0.786 

0.733 

0.673 

0.608 

0.538 

0.464 

0.389 

0.315 

0.247 
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0.135 
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Figure  26.  Air-sea  Interface  transmittance  as  a function  of  sun  or  moon  zenith  angle 
and  surface  wind  speed,  v. 
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For  diffuse  or  uniform  radiation  incident  on  the  sea  surface,  we  use  the  approxi- 
mation in  section  2.3.4  (which  neglects  wave  effects)  and  take  r^w  j = 0.83. 


Tawls  = 0'83’  ropt>l°; 

Tawlm  = O-83’  Topt  ^ *0’ 
and  for  blue  sky, 

raw  IB  = 0.83,  all  values  of  rQpt, 

while  for  stellar/zodiacal  light, 

Taw  1 Z = ^ 83>  values  Topt 
3.3.5  Air-Water  Interface  Angular  Effects 


(N-9) 


(N-10) 


(N-l  1) 


(N-12) 


The  wave  slopes  on  the  sea  surface  cause  an  overall  increase  in  the  beam  divergence 
of  an  incident  beam,  or  equivalently,  the  apparent  angular  size  of  the  source  as  viewed  from 
an  underwater  point  of  view.  With  regard  to  the  background  sources,  only  the  sun  and 
moon  for  clear  weather  conditions  (rQpt  < 10)  will  be  appreciably  affected. 

Again,  using  the  Karp  model  discussed  in  section  2.3.5, 


^sujnu  « 0.0l03v%(rop|<l0); 


(N-13a) 


for 


^!uw 


For  all  r 


for 


opt’ 

A flB.Z 
aw 


A*awD 

^aw^ 


= rms  induced  half-angle  spread  for  the  sun,  and 

= rms  induced  half-angle  spread  for  the  moon. 

= surface  wind  speed  in  knots. 

= 0 ; 

= effect  on  blue  sky  source,  and 
- effect  on  stellar/zodiacal  source. 


(N-14) 


Also 


mu  = o,  (r 


opt ’ 


10), 


since  the  light  is  diffuse  after  emerging  from  the  thick  clouds. 


(N-13b) 
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Table  1 7 and  figure  27  evaluate  equation  (N-I3a)  for  V in  knots  (and  metres  per 
second). 


Table  17.  Rms  air-water  interfae  induced  half-angle  effects 
(r0pt  < 10)  sun  a 'i'  i moon. 


v,  Wind  Speed 

Aflaw  su  or  mu 

Knots 

Metres/Sec 

Milliradians 

Degrees 

0 

0 

0 

0 

2 

1.03 

14.6 

0.84 

4 

2.06 

20.7 

1.18 

B 

4.12 

29.2 

1.67 

14 

7.21 

38.6 

2.21 

20 

10.3 

46.2 

2.65 

26 

13.4 

52.6 

3.0 

32 

16.5 

58.4 

3.35 

38 

19.6 

63.6 

3.64 

SURFACE  WIND  SPEED.  matr»>/iec 


SURFACE  WIND  SPEED  (know),  V 


Figure  27.  Rms  air-water  interface  effect  as  a function  of  wind  speed  v. 

Since  the  full  angular  subtense  of  the  sun  (and  the  moon)  is  *0.5  degrees,  this  effect 
will  substantially  increase  its  apparent  size.  The  relative  contribution  of  A0aw  to  the  distri- 
bution of  noise  radiance  at  the  receiver  is  discussed  in  section  3.3.8.  Except  for  the  clearest 
water  it  is  a small  effect,  and  so  the  impact  of  neglecting  zenith  angle  effects  and  dissimilar 
wave  slopes  in  the  downward  and  crosswind  direction  may  be  negligible.  We  therefore 
adopt  this  model  until  better  information  is  available. 
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DEGREES 


3.3.6  Relative  Surface  Foam  Coverage 


The  energy  transmission  of  the  air-water  interface  is  composed  of  two  factors: 
raw  = (t^w1)  X (raw2), 
for 

r’,w  total  e nervy  transmission  of  the  air-water  interface 

rawl  “ iir-water  interface  transmission  due  to  index  of  refraction  discontinuity; 

and, 

Taw2  = air'water  interface  transmission  due  to  foam  and  streaks  on  the  water  surface. 

This  section  treats  r^.T.  while  r'awj  has  been  treated  in  section  3.3.4. 

The  surface  foam  coverage  and  its  effects  are  taken  to  be  independent  of  the  noise 
source  and  cloud  conditions.  As  discussed  in  section  2.3.6,  for  a foam  albedo  = 1 , 

r'aw2=  1 -0-2  (1(T5))  v3  J,  v ^9  m/sec  , 
and 

Taw2=  1 - ( 1 -2  < 10~5))  v 3 3 (0.225v -0.99),  v > 9 m/sec 
for 

v = surface  wind  speed  in  metres/sec. 

Equation  (N-l  5a,  b)  is  evaluated  in  table  18  and  figure  28  for  v in  knots  (and  metres/second). 

Although  this  model  neglects  zenith  angle  effects  we  shall  adopt  it  pending  further 
experimental  work. 


(N-l  5a) 

(N-l  5b) 


Table  18.  Air-water  energy  transmission  due  to  surface  foam  and  streaks 
(assuming  a foam  streak  albedo  = 1 ). 


v,  Wind  Speed 

raw2 

knots 

met  res/  sec 

0 

0 

1 

2 

1.03 

I 

4 

2.06 

l 

8 

4.12 

1 

14 

7.21 

0.99 

20 

10.3 

0.96 

26 

13.4 

0.87 

32 

16.5 

0.66 

38 

19.6 

0.25 
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SURFACE  WIND  SPEED  (knott),  v 


Figure  28.  Foam/streak  surface  coverage  transmission  vt  .sus  surface  v ind  speed. 

3.3.7  Water  Energy  Transmission 

The  energy  transmission  of  the  water  is  denoted  by  r'w.  The  angularly  localized 
noise  sources  (sun  and  moon)  behave  similarly  to  the  signal  energy  transmission  discussed 
in  section  2.3.7,  thus  we  take: 


J 

I <kiDi> 

r = exn  - 

i=l 

'wsu  ^ 

cos  ^sWu 

for 

<u  =sin'’ 

1 Sin  4>su 

(N-16) 


I D‘-D; 

i=l 


where, 


= diffuse  attenuation  coefficient  for  the  i’th  water  layer, 
= thickness  of  the  i’th  water  layer, 

= eceiver  depth, 

= sea-water  index  of  refraction, 
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<p*u  = in-water  solar  zenith  angle, 
0SU  = in-air  solar  zenith  angle; 


and, 


i=  1 


(N-17) 


Here, 

^mu  = in*v  lunar  zenith  angle; 

0mu  = - ,iar  zenith  angle. 

For  the  diffuse  or  uniform,  background  sources,  such  as  the  blue  sky  and  the  stars, 
we  assume  there  is  no  zenith  angle  dependence;  and  so. 


j 

r'wb  = exp  - (kjDj) 
1=1 


and 


T WZ  = CXP  ’ 


(N-18) 


(N-19) 


We  do  .101  distinguish  here  between  clear  sky,  thin  cloud,  and  thick  cloud  models. 
This  model  is  uncertain,  in 

(1)  the  values  of  kj  to  use, 

(2)  the  values  of  Dj,  and 

(3)  its  applicability  in  very  clear  water  and/or  at  shallow  receiver  depths. 

It  is  the  best  model  available  now  and  it  will  be  revised  when  better  information 
becomes  available. 


-~1 


3.3.8  Water  Distribution  of  Radiance 

There  is  no  experimentally  verified  expression  for  the  in-water  distribution  of  back- 
ground radiance  as  a function  of  source  character,  source  zenith  angle,  water  properties 
and  receiver  depth.  As  discussed  in  section  2.3.8,  we  therefore  adopt  the  expression 

i - ( »w-\2 , 

\ s,n  <t>0  ) 

as  an  estimate  of  the  angular  distribution,  with 


0W  - in-water  angle  measured  from  the  axis,  or  principal  ray  of  the  noise  source. 
Here,  0g  is  related  to  the  half-power  point  of  the  received  radiance  by  the  equation, 


1 - (cos  01^)  - i 

j 3 sin2  <Pq 

cos  (0i/2)  sin2  (0y2)  + 2 cos  (014)  - 2 

2 1 

1 - cos  <p0 

3 sin2  0q 

cos  0Q  sin“  00  + 2 cos  0g  - 2 

(N-20a) 


Equation  (N-20a)  is  evaluated  in  table  19.  Values  between  those  shown  are  obtained  by 
linear  interpolation. 


Table  19.  Relation  of  radiance  zero  point,  0q,  and 
received  radiance  half-power  point,  for 
1 - (sin  0w/sin  0g)2  radiance  distribution. 


0(4  (degrees) 

01/4  (degrees) 

3.8 

5 

7.6 

10 

11.4 

15 

15.2 

20 

19.0 

25 

22.7 

30 

26.5 

35 

30.2 

40 

33.9 

45 

37.5 

50 

41.1 

55 

44.6 

i 

60 

48.1 

65 

51.6 

70 

54.9 

75 

58.2 

80 
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Again,  assuming  that  the  in-air  incident  beam  spread,  air-water  beam  spread,  and  in- 
water scattering  induced  spread  are  statistically  independent  effects,  we  adopt  the  NOSC 
(ref  7)  model: 


~~  V*  + *aw  + ^ ’ 


for  all  four  out-of-water  background  sources.  For  solar  and  lunar  sources, 


and 


(N-20b) 


(N-21a) 


*W  *Si a!1  TOpt  ; 

cos 

while,  for  the  distributed  background  sources  of  blue  sky  and  stellar/zodiacal  light, 


(N-22a) 


fw  = sD  , all  r 


opt 


(N-23a) 


Here, 


* 

-I 


fljj  * mean  square  single  scattering  angle  in  water, 

s * scattering  coefficient  in  water, 

D = receiver  depth, 

^su  = in-water  solar  zenith  angle,  and 

^mu  = in-water  lunar  zenith  angle. 

Again,  for  solar  and  lunar  sources, 

n 

(N-2 1 b),  (N-22b) 
(N-21c),  (N-22c) 

V = surface  wind  speed  (knots),  as  discussed  in  section  4.3.5  For  the  distributed 


sources, 


faw  = 0 , all  r. 


opt 


(N-23b) 
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Finally,  for  the  sun  and  moon, 


-] 


(N-2ld) 
(N-2 1 e) 


/; 


and 


> T opt  ^ 

> 10  . 


Here, 


n = water  index  of  refraction, 

d$j2  = half  the  angular  substense  of  the  sun  (~( 1/4)°),  and 
Bmj2  ~ half  the  angular  substense  of  the  moon  (~l/4)°). 


(N-22d) 

(N-22e) 


These  equations  have  been  discussed  and  derived  in  section  2.3.8.  For  the  distributed 


sources. 


f=(33.8°)2  , all  r 


opt 


In  general,  the  receiver  will  be  directly  viewing  the  signal,  while  the  background  source 
enters  at  an  off-axis  angle.  Then,  the  fraction  of  the  noise  radiance  which  enters  the  receiver 
is  given  by 


for 


0W  = cos“*  jeos  0W  cos  5 + sin  0W  sin  6 sin  6 j ; 

where, 

- half-angle  of  the  receiver  field  of  view, 

5 = off-set  angle  between  axis  of  noise  source  and  receiver  optical  axis. 
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This  expression  will  be  used  further  in  section  3.2.10  and  3.2.1 1. 

Equation  (N-24a)  applies  to  the  background  sources  in  thin  cloud  conditions.  Under 
thick  cloud  conditions,  both  signal  and  background  will  appear  to  arrive  from  the  zenith, 
and  so  5 -*  0.  For  this  case, 


f(*o-  ~ 


cos  0R  sin^  0R  + 2 cos  - 2 

t -cosoR- 

3 stnz  <(>q 

cos  0q  sir/  + 2 cos  - 2 

I - --  " ■ 

3 sinz  0q 

(N-24b) 


3.3.9  Detection  Bandwidth 


The  required  electrical  detection  bandwidth  to  optimally  detect  the  pulses  discussed 
in  section  2.3.9  is  not  known  at  present.  In  lieu  of  such  a result  we  assume: 

1)  The  receiver  has  foreknowledge  of  the  expected  pulse  width,  and 


2) 


B = - 


0.4 


(2.45  tm) 


(N-25) 


for, 

B = electrical  detection  bandwidth, 

tm  = time  at  which  pulse  peak  value  occurs  after  pulse  start,  for  a pulse  shape 
f(t)  * t exp  - (t/tm). 

For  gaussian  shaped  pulses  and  the  detection  filter,  equation  (N-25)  is  the  nearly 
optimum  match.*  As  further  work  is  done  in  the  area  of  the  real  pulses  to  be  expected 
here,  (N-25)  may  be  revised. 

Equation  (N-25)  is  evaluated  in  table  20  and  figure  29  for  the  pulse  widths  and 
optical  thicknesses  developed  in  section  2.3.9. 


*H.P.  Weslman,  Editor,  Reference  Data  for  Radio  Engineer,  Fifth  Edition,  p.  29-5,  (H.W.  Sams  A Co., 
New  York,  1969). 
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3.3.10  Average  Background  Power  Due  to  Sunlight 

The  average  optical  background  power  in  the  receiver  due  to  the  sun  is  given  by  an 
equation  analogous  to  that  developed  in  section  2.3. 10  for  the  received  optical  signal 
energy.  We  therefore  take 

Psu  = (spectral  irradiance  above  the  atmosphere  due  to  sun), 

X (clear  atmospheric  transmission)  X (cloud  transmission) 

X (cloud  to  water  transmission) 

X (air-water  interface  transmission)  X (water  transmission) 

X (receiver  optics  transmission) 

X (fraction  of  incident  radiance  within  receiver  field-of-view) 

X (area  of  the  receiver).  (N-26) 


We  also  take 


sX 


Dopt 

T'R 

d 


(rrd“/4) 


' cw 


'aw 


w 


exo-atmospheric  spectral  solar  irradiance, 
spectral  bandwidth  of  receiver, 
receiver  optics  transmission, 
diameter  of  receiver  aperture, 

area  of  receiver,  (N-27) 

clear  atmosphere  energy  transmission  as  discussed  in  section  2.3.1, 

cloud  energy  transmission  as  discussed  in  section  2.3.2, 

cloud  to  water  energy  transmission  as  discussed  in  section  2. 3. 3, 

air-water  interface  energy  transmission  as  discussed  in  sections  2.3.4 
and  2.3.6, 

water  energy  transmission  as  discussed  in  section  2.3.7,  and 


1(0*)  = water  radiance  distribution,  as  discussed  in  section  2.3.8. 


Therefore,  the  received  optical  power  is  given  by 


Psu  * HsX  Bopt  ^d2/4>  7R>  Ta  Tc  rcw  raw  Tw  f(*o-  <N’28) 

The  fraction  of  the  incident  radiance  within  the  receiver  field-of-view  is  given  by  (for  perfect 
alignment  between  beam  axis  and  receiver  axis): 


2n 

Jo  i 

f/O'i  - u 

0R 

f I(0W)  d Cl 
° a,,  < a. 

f(n)  2rr 

Jo  J 

4>o  ' K 0 

f K0w)dn 

0 

4 


|| 


J 

i 

i 

i 

.1 

i 

i 

t 

1 

i 

! 
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for 


= half-angle  of  the  receiver  field-of-view 
0O  = off-axis  angle  at  which  incoming  radiance  equals  zero. 
Using  the  model  adpted  in  section  2.3.8, 


(N-30) 


(N-31 ) 


3.3.1 1 Average  Background  Power  Due  to  Moonlight 

The  average  optical  background  power  in  the  receiver  due  to  the  moon  is  completely 
analogous  to  that  for  the  sun  discussed  in  section  3.3. 10.  We  therefore  take 

Pmu  = (spectral  lunar  irradiance  outside  the  atmosphere) 

X (area  of  the  receiver) 

X (clear  atmospheric  transmission)  X (water  transmission) 

X (cloud  to  water  transmission) 

X (air-water  interface  transmission)  X (water  transmission) 

X (receiver  optics  transmission) 

X (fraction  of  incident  radiance  within  receiver  field-of-view).  (N-32) 

We  take 

Hmux  3 exo-atmospheric  lunar  irradiance, 

Bopt  = spectral  bandwidth  of  receiver, 

■y^  = receiver  optics  transmission, 

d ® diameter  of  receiver  aperture, 

ra  = clear  atmosphere  energy  transmission  as  discussed  in  section  2.3. 1 , 

rc  * cloud  to  energy  transmission  as  discussed  in  section  2.3.2, 

rcw  = cloud  to  water  energy  transmission  as  discussed  in  section  2.3.3, 

= air-water  interface  energy  transmission  as  discussed  in  sections  2.3.4 
aW  and  2.3.6, 

rw  - water  energy  transmission  as  discussed  in  section  2.3.7,  and 
1(0W)  - water  radiance  distribution  as  discussed  in  section  2.3.8. 


Therefore,  the  received  optical  energy  is  given  by 


^mu  = ^muX  ®opt  ('  4 ) ra  Tc  Tcw  Taw  Tw  f^o’ 


being  the  area  of  the  detector: 


2t  6 d 

f f Htw) d a 

f(«)  = 

f t°  K#w)dn 

' o ' o 


aR  ^o 


0R  = half-angle  of  the  receiver  field-of-view, 

<t>Q  = off-axis  angle  at  which  incoming  radiance  equals  zero. 

Using  the  model  adopted  in  section  2.3.8, 


f(fl)  = f(^w,  0R) 


for  0 R < <t>0  ■ 


f(0Q.  0r)  = 1 , for0R>tfo  ; 

3.3.12  Average  Background  Power  Due  to  Blue  Skylight 


(N-33) 


(N-34) 


(N-35) 


1 -COS0R  ~(<t>0)~2 

20R  sin  0r-(0r-2)  cos  0R  - 2 

1 - cos  <t>0  - ( <t>0)~ 2 

2<t>l  sin  <j>Q  - ( <t>l  - 2)  cos  <t>0  - 2 

(N-37) 


The  average  optical  background  power  in  the  receiver  due  to  the  blue  skylight  is 
partially  analogous  to  that  for  the  sun  and  moon  discussed  in  sections  3.3. 10  and  3.3. 1 I. 

We  therefore  take,  for  the  average  optical  background  power  due  to  blue  skylight: 

PBS  = (spectral  radiance  at  receiver  aperture  due  to  the  blue  sky) 

X (teceiver  optics  transmission)  X (receiver  area) 

X (receiver  optical  filter  bandpass)  X (receiver  solid  angle) 

X (fraction  of  incident  radiance  within  the  receiver  field  of  view).  (N-38) 
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We  also  take 


- receiver  optics  transmission, 
* receiver  aperture  diameter. 


= area  of  receiver  aperture, 

= receiver  optical  filter  bandpass, 

= half-angle  of  the  receiver  field-of-view, 


2jr(l  - cos  0g)  = receiver  field-of-view  solid  angle, 

Lg§  = spectral  radiance  at  receiver  aperture  due  to  the  blue  skylight, 

C (<t>0,  5)  = fraction  of  incident  radiance  within  receiver  field  of  view,  and 


PBS"LBS  (7R  (-4-)  ) Bopt  2jr(*  ~cos0r)  1v(0o,  0r,6>. 


(N-39) 


Because  the  blue  skylight  is  a cw  source,  we  use  the  energy  transmission  formalism 
to  develop  the  expression: 

Lbs  = (clear  sky  exo-atmospheric  effective  radiance) 

X (clear  atmospheric  transmission)  X (cloud  energy  transmission) 

X (cloud  to  water  energy  transmission) 

X (air-water  interface  energy  transmission)  X (w  ater  energy  transmission); 


and  we  use 


(N-40) 


clear  sky  exo-atmospheric  effective  radiance. 


T3  = clear  atmospheric  transmission  as  discussed  in  section  3.3.1, 

TCB  = cloud  energy  transmission  as  discussed  in  section  3.3.2, 

Tcw  = cl°ud  (o  water  energy  transmission  as  discussed  in  section  3.3.3, 

TawB  = a'r'water  energy  transmission  as  discussed  in  sections  3.3.4  and  3.3.6,  and 
rwB  = water  energy  transmission  as  discussed  in  section  3.3.7. 

Gathering  the  expressions,  we  find 


^BS  Ta  rcb  rcw  rawB  TwB 


(N-41) 
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Again,  the  fraction  of  blue-sky  radiance  within  the  receiver  field  of  view  is  given  by 


u 
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o n 

f dO  f d0w  sin 

neR,»0«>* — ^ 


/ « f 

J 0 Jo 


ft>  0 


d0w  sin  <t> 


for 


(N-42) 


(N-43) 


and 

6 = off-zenith  pointing  angle  of  the  receiver  axis, 

while 

4>q  = off-zenith  angle  at  which  the  blue  sky  ndiance  goes  to  zero. 

These  expressions  will  be  used  in  the  development  of  the  noise  equivalent  optical 
power  expression  in  section  3.3.15. 

3.3.13  Average  Background  Power  Due  to  Stellar/Zodiacal  Light 

The  average  optical  background  power  in  the  receiver  due  to  the  nighttime  distri- 
buted sources  of  stellar  or  zodiacal  light  follows  the  patterns  established  in  the  previous 
three  sections.  We  take,  for  the  average  optical  background  power  due  to  these  sources: 

Pj  = (spectral  radiance  at  receiver  aperture  due  to  the  stellar/zodiacal  light) 

X (receiver  optics  transmission)  X (receiver  area) 

X (receiver  optical  filter  bandpass)  X (receiver  sold  angle) 

X (fraction  of  incident  radiance  within  the  receiver  field-of-view).  (N-44) 

receiver  optics  transmission, 
receiver  aperture  diameter, 

area  of  receiver  aperture, 

receiver  optical  filter  bandpass, 
half-angle  of  the  receiver  field-of-view, 
receiver  field-of-view  solid  angle, 

spectral  radiance  at  receiver  aperture  due  to  the  stcller/zodiacal 
light, 


We  again  take 

d 

ffd* 

4 

®opt 

0R 

2jt(!  - cos  0 r)  = 
LZS 
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r(t0’  0r> 

*o 


Then, 


* fraction  of  incident  radiance  within  the  receiver  field-of-view,  and 

= off-axis  angle  at  which  the  received  radiance  goes  to  zero,  as 
discussed  in  section  3.3.8. 


pz = Lzs 


(2rr(l  -cos0R))  f(0o,  eR<  5> 


(N-45) 


Because  this  background  source  is  cw,  we  use  the  energy  transmission  formalism  to 
develop  the  expression: 

Lzs  * (stellar/zodiacal  light  clear  sky  effective  exo-atmospheric  radiance) 

X (clear  atmospheric  transmission)  X (cloud  energy  transmission) 

X (cloud  to  water  energy  transmission) 

X (air-water  interface  energy  transmission) 

X (water  energy  transmission).  (N-46) 

We  also  use, 

L2  = stellar/zodiacal  light  clear  sky  effective  exo-atmospheric  radiance, 
r'a  = clear  atmospheric  transmission  as  discussed  in  section  3.3. 1 . 

r'cz  = cloud  energy  transmission  as  discussed  in  section  3.3.2, 

t'cw  = cloud  to  water  energy  transmission  as  discussed  in  section  3.3.3, 

Tawz  = air-water  interface  energy  transmission  as  discussed  in  sections  3.3  4 
and  3.3.6,  and 

Twz  * water  energy  transmission  as  discussed  in  se^ion  3.3.7. 


Gathering  the  expressions  we  find: 


^ZS  ~ ®"Z  Ta  rcz  Tcw  rawz  rwz  • (N-47) 

Again,  the  fraction  of  stellar/zodiacal  radiance  within  the  receiver  field-of-view  is 
given  by 
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(N-48) 


for 


0W  = cos  - ! J wcos  0W  cos  5 + sin  0W  + sin  6 sin  0 j 
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and 


£ 

j- 


i 


& = Off-zenith  pointing  angle  of  the  receiver  axis. 

These  expressions  will  be  used  in  the  development  of  the  noise  equivalent  optical 
power  expression  in  section  3.3.15. 


3.3.14  Average  Background  Power  Due  to  Bioluminescence 


The  final  source  of  optical  background  power  is  the  local  bioluminescent  sources 
which  are  stimulated  to  emit  by  the  submarine  motion,  or  other  disturbances  in  the  water 
This  is  modeled  in  a slightly  different  way  than  the  previous  four  sources,  and  cloud  and 
water  properties  have  only  an  indirect  effect  on  this  source  strength.  We  therefore  write 
for  the  average  background  power  due  to  bioluminescence. 


^BL  = (sPectra*  irradiance  at  receiver  aperture  due  to  bioluminescence) 

X (receiver  optics  transmission)  X (receiver  area) 

X (receiver  optical  filter  bandpass); 

and  we  set 

= spectral  iiradiance  at  receiver  aperture  due  to  bioluminescence, 
= receiver  optics  transmission, 

= diameter  of  receiver  aperture. 


(N-50) 


LBL 

*R 

d 


7rd 

4 


= areu  of  receiver  aperture,  and 
Bopt  = optical  filter  bandpass. 


Therefore, 

PBL  = LBL  (fR 


(x)) 


opt 


(N-51) 


3.3. 1 5 Noise  Equivalent  Optical  Power  Dependence  on  Noise  Sources 


In  general,  a direct  detection  optical  communication  system  has  four  independent 
noise  contributions  which  include  thermal  (or  amplifier)  noise,  dark  current  detector  noise, 
signal  shot  noise,  and  background  shot  noise.  These  are  noise  sources  insofar  as  they  gener- 
ate fluctuations  in  the  electrical  current  present  in  the  detection  system.  It  is  conventional 
to  write  the  “noise”  as  the  1 -a  point  of  the  fluctuating  electrical  current,  assuming  the 
noise  sources  add  independently  and  are  steady  in  character. 

Because  we  have  derived  a signal  level  in  terms  of  the  instantaneous  received  optical 
power,  it  is  appropriate  to  describe  the  noise  components  in  terms  of  a noise  equivalent 
(optical)  power,  as  derived  from  the  post-detection  electrical  power. 


Ultimately,  we  write,  for  a photomultiplier  tube  type  of  detector. 
NtPto.=  [NEPJ^  + NEp*jc  + NEPjs  + NEPg  J * ; 


(N-52) 
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NEPtot 
NEPth 
NEP  dc 
NEPSS 

NEPg 


total  noise  equivalent  (optical)  power  due  to  all  sources, 


noise  equivalent  optical  power  due  to  thermal  or  amplifier  noise, 

noise  equivalent  optical  power  due  to  photo-detector  dark  current, 

noise  equivalent  optical  power  due  to  shot-noise  generated  by  the  signal,  and 

noise  equivalent  optical  power  due  to  shot-noise  generated  by  the 
background. 


1 


Then. 

(N-53) 

lor, 

(kT)  = thermal  noise  energy  = (Boltzman’s  constant)  X (absolute  temperature), 

B = electrical  detection  bandwidth  as  discussed  in  section  3.3.9, 

G = detection  gain, 

9 = photo  surface  quantum  efficiency, 

e = charge  on  the  electron, 

]u>  = energy  per  signal  photon, 

r?e/hv  = photo-surface  responsivity, 

= load  resistance,  and 

Fa  = amplifier  noise  figure. 

For  the  dark  current  contribution, 

(N-54) 

for, 
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For  the  signal  shot  noise  contribution, 


2 e B F (n<*/h»>)  G^PR  Rl 
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lA 

NEPSS  = 

2 e B FPr 

r}e/hv)2  R|_ 

(tje/hy) 

for  PR  = peak  received  optical  signal  power  at  the  photo  surface,  as  discussed  in 
section  2.3.10. 


Finally,  the  cw  background  contributes 


nepb  = 

'2eBFG2(tjeW  RL|S  P^j 
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2 e B F Z P}} 
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s PB  " psu  + Pmu  + PBS  + PZ  + PBL 
i 


(N-55) 


(N-56) 


(N-57) 


where  (see  fig.  20,  page  47/48) 

Psu  = average  background  power  due 
Pmu  = avera8e  background  power  due 
Pbs  = average  background  power  due 

- average  background  power  due 
section  3.3.13,  (N-45); 


to  sunlight  as  discussed  in  section  3.3. 10,  (N-27); 
to  moonlight  as  discussed  in  section  3.3. 1 1 , (N-33); 
to  blue  skylight  as  discussed  in  section  3.3. 1 2,  (N— 39) ; 
to  stcllar/zodiacal  light  as  discussed  in 


and 

pBL  = average  background  power  due  to  biolumincscence  as  discussed  in 
section  3.3. 14,  (N-5  1 ). 


Two  comments  are  in  order  at  this  point. 


1.  If  the  signal  shot  noise  dominates  the  noise  components,  the  formulation  should 
be  re-examined  to  insure  that  enough  photo-electrons  are  being  generated  to  make 
it  applicable. 

2.  Not  all  the  average  background  contributors  will  be  present  at  any  one  time, 
which  will  be  accounted  for  in  the  time-of-day  modeling  of  the  respective 
spectral  radiances. 
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GLOSSARY 


This  glossary  defines  all  the  English  and  Greek  symbols  used  in  the  Naval  blue-green 
downlink  model.  The  English  symbols  are  listed  alphabetically,  followed  by  the  Greek 
symbols. 

ENGLISH  SYMBOLS 

A Ag  = Energy  to  instantaneous  power  normalization  parameter 

B b = Effective  clear  atmosphere  optical  thickness 

B = Electrical  detection  bandwidth 
Bopt  = Receiver  optical  filter  bandpass 

C Cf  = Fraction  of  sea-surface  covered  by  foam  and  streaks 
c = Speed  of  light 

D D = Receiver  depth 

Dj  Thickness  of  i’th  water  layer 

d = Diameter  of  receiver  aperture 

E = Total  received  energy  per  pulse 

exp  = Exponential 

Ep  = Transmitted  energy  per  pulse 

e = Charge  on  the  electron 
E2  = Exponential  integral 

F fw  - Water  contribution  to  received  beam  half-angle 

fcw  = Air-water  interface  contribution  to  received  beam  half-angle 
f = Atmospheric  contribution  to  received  beam  radiance 

a 

fwj  = Contribution  of  i’th  water  layer  to  received  beam  half  angle 
fit)  = Received  pulse  shape 

F (4>  o-  fR)  = Fraction  of  incident  radiance  within  receiver  field  of  view 
F(0O.  #R,  $su)  = Fraction  of  incident  radiance  within  receiver  field  of  view 

F = Excess  noise  in  detector  gain 
Fa  = Amplifier  noise  figure 

G G = Detection  gain 

H H = Distance  from  cloud  base  to  water  surface 

h^  = Energy  per  signal  photon 
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I I(0W)  = Water  radiance  distribution 

* Dark  current  at  the  photo-cathode 

J j = Number  of  water  layers  present  from  surface  to  submarine  receiver 

K k = Diffuse  attenuation  coefficient  of  the  water 

kj  = Diffuse  attenuation  coefficient  of  i’th  water  layer 

(kt)  = Thermal  noise  energy  = (Boltzmann’s  constant)  X (absolute  temperature) 

L Lgjj  = Spectral  radiance  at  receiver  aperture  due  to  the  sun 
L§  = Exo-atmospheric  effective  solar  radiance 
Lmu  “ Spectral  radiance  at  receiver  aperture  due  to  the  moon 
Lm  = Exo-atmospheric  effective  lunar  radiance 
Lfis  “ Spectral  radiance  at  receiver  aperture  due  to  blue  skylight 
Lr  = Clear  sky  exo-atmospheric  effective  radiance  due  to  blue  skylight 
Lzs  = Spectral  radiance  at  receiver  aperture  due  to  stellar  and  zodiacal  light 

= Clear  sky  exo-atmospheric  effective  radiance  due  to  stellar  and  zodiacal  light 
Lql  = Spectral  irradiance  at  receiver  aperture  due  to  bioluminescence 

N n * Water  index  of  refraction 

NEPtot  = Total  noise  equivalent  optical  power  due  to  all  sources 

nepTh  = Noise  equivalent  optical  oower  due  to  thermal  or  amplifier  noise 
NEPdc  = Noise  equivalent  optical  power  due  to  photo-detector  dark  current 
NEP<;s  = Noise  equivalent  optical  power  due  to  shot-noise  generated  by  the  signal 
nepb  = Noise  equivalent  optical  power  due  to  shot-noise  generated  by  the  background 

P Pr(*)  = Instantaneous  received  signal  power 

Psu  = Average  optical  power  at  receiver  due  to  the  sun 
Pmu  = Average  optical  power  at  receiver  due  to  the  moon 
PBS  = Average  optical  power  at  receiver  due  to  the  blue  sky 
Pz  - Average  optica)  power  at  receiver  due  to  stellar  and  zodiacal  light 
PgL  = Average  optical  power  at  receiver  due  to  bioiuminescence 
= Peak  received  signal  power 

Q q = Parameter  describing  ability  of  satellite  transmitter  to  correct  for  zenith  angle 
spot  spreading.  0 < q < 1 
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R R = Range  from  satellite  to  submarine 
R^  = Load  resistance 

5 s = Water  scattering  coefficient 

Sj  = Scattering  coefficient  of  i’th  water  layer 

T T = Geometrical  thickness  of  the  cloud 
t - Time 

tm  - Time  after  pulse  start  at  which  peak  value  occurs 
Atw  = Pulse  width  due  to  water  portion  of  the  path 
Atcw  = Pulse  width  due  to  cloud  to  water  portion  of  the  path 
Atc  = Pulse  width  due  to  cloud  portion  of  the  path 

U 

V v = Surface  wind  speed 
W 

X 

Y 
Z 

GREEK  SYMBOLS 

The  Greek  symbols  are  listed  in  order  according  to  the  Greek  alphabet:  a,  0,7,  5, 
e,  f,  V’  P.  o,  i r,  p,  a,  r,  v,  0,  x,0.  u>- 

a 

0 

7 7j  = Transmitter  optics  transmission 
7^  = Receiver  optics  transmission 

6 6 = Off-set  angle  between  receiver  optical  axis  and  axis  of  the  incoming  light 
e 

f 

n 
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6 * Cloud  particle  mean  scattering  angle 
<cos  0>  = Mean  cosine  of  the  in-cloud  scattering  angle 
0-p  = Full  angle  exp  (-2)  transmitter  beamwidth 

AS  aw  = ^ms  half-angle  additional  signal  beam  divergence  due  to  wave  action 

2 

= Mean  square  single  scattering  angle  in  water 

2 

0gij  = Mean  square  single  scattering  angle  in  water  for  the  j’th  water  layer 
= Half  angle  of  the  receiver  field  of  view 

AS  aw  = Rms  half-angle  air-water  interface  induced  spread  for  the  sunlight 

ASaw  = Rms  half-angle  air-water  interface  induced  spread  for  the  moonlight 
S$/l  = Half  the  angle  subtended  by  the  sun 
Sm/l  = Half  the  angle  subtended  by  the  moon 

AS®w  = Rms  half-angle  air-water  interface  induced  spread  for  the  blue  sky  light 

7 

AS^w  = Rms  half-angle  air-water  interface  induced  spread  for  the  stellar/zodiacal  light 


<7C  = Mean  extinction  coefficient  of  the  cloud 

Ta  = Signal  clear  atmospheric  energy  transmission 

tc  - Signal  cloud  energy  transmission 

Topt  = Ophcal  thickness  of  the  cloud 

rcw  = Signal  cloud  to  water  energy  transmission 

raw  = Signal  total  energy  transmission  of  air-water  interface 
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rawl  = Signal  air-water  interface  energy  transmission  due  to  the  index  of  refraction 
discontinuity 

Taw2  = Signal  air-water  interface  energy  transmission  due  to  foam  and  streaks  on 
the  sea  surface 

rw  = Signal  water  energy  transmission 

r'a  = Background  clear  atmospheric  energy  transmission 

Tcb  = Cloud  energy  transmission  of  the  blue  skylight 
= Cloud  energy  transmission  of  direct  sunlight 
rcm  = Cloud  energy  transmission  of  the  direct  moonlight 
t'cz  = Cloud  energy  transmission  of  the  stellar  and  zodiacal  light 
t'cv/  = Background  cloud  to  water  energy  transmission 

raw  = Total  background  energy  transmission  of  the  air-water  interface 

Tawl  = Background  air-water  interface  transmission  due  to  index  of  refraction 
discontinuity 

raw2  = Background  air-water  interface  transmission  due  to  foam  and  streaks  on 
the  sea  surface 

Tawis  = Solar  air-water  interface  transmission  due  to  index  of  refraction  discontinuity 

Tawim  = Lunar  air-water  interface  transmission  due  to  index  of  refraction  discontinuity 

Tawib  = Blue-skylight  air-water  interface  transmission  due  to  index  of  refraction 
discontinuity 

Tawiz  = Stellar  and  zodiacal  light  air-water  interface  transmission  due  to  index  of 
refraction  discontinuity 

rwsu  = Solar  water  energy  transmission 
rwmu  = Lunar  water  energy  transmission 
= Blue  sky  water  energy  transmission 
V Twz  = Stellar  and  zodiacal  light  water  energy  transmission 

<t>  0S  = Signal  in-air  zenith  angle 

0^  = Signal  in-water  zenith  angle 

00  = Off-axis  angle  at  which  in-water  radiance  goes  to  zero 

0i/2  = Half-power  angle  of  the  received  signal  beam  radiance 

0SU  = Solar  in-air  zenith  angle 

0mu  = Lunar  in-air  zenith  angle 

N01/ 2 = Half-power  angle  of  the  background  radiance 
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